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FREE  MOLECULE  FLOW  THEORY  AND  ITS  APPLICATION 
TO  THE  DETERMINATION  OF  AERODYNAMIC  FORCES 

ABSTRACT 

The  purpose  of  this  report  is  twofold.  It  is  designed  as  an  introduction 
to  free  molecule  flow  theory  for  those  with  no  background  in  rarefied  gas 
dynamics  or  kinetic  theory,  and  it  presents  an  exact  method  for  the  determi¬ 
nation  of  the  forces  and  moments  on  bodies  in  a  free  molecule  flow.  Approx¬ 
imate  methods  are  evaluated. 

Comoarison  with  the  exact  method  of  this  report  shows  that  for  molecu¬ 
lar  speed  ratios  at  satellite  altitudes  the  approximate  theories,  which  do  not 
completely  account  for  the  random  motion  of  the  impinging  molecules,  are 
satisfactory  for  the  prediction  of  C^  and  C  For  elongated  bodies  at 
small  angles  of  attack,  the  approximate  methods  are  shown  to  be  inadequate 
for  the  prediction  of  C^.  Tables  of  certain  functions  requisite  to  the  detei  - 
mination  of  free  molecule  aerodynamic  coefficients  are  presented. 
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i,  j,  h 

Kn 
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-  area 

-  average  speed  of  a  molecule 

-  axial  force  coefficient 

-  drag  coefficient 

-  lift  coefficient 

-  pitching  moment  coefficient 

-  normal  force  coefficient 

-  most  probable  random  speed  of  the 
moleculesj  -  J Z  R  T' 

-  Free  molecule  flow  functions  for  a  flat 
plate  with  both  sides  exposed  to  the  flow 

-  P'ree  molecule  flow  functions  for  a  flat 
plate  with  one  side  exposed  to  the  flow 

-  F  ree  molecule  flow  functions  for  a 
cylinde  r 

1^®  .  F  ree  molecule  flow  functions  for  a  cone 


-  Free  molecule  flow  functions  for  a 
spherical  segment 

-  unit  vectors 

-  Knudsen  number 

-  direction  cosine  between  the  local  x  axis 
and  the  direction  in  which  the  force  is 
desired 


1  -  direction  cosine  between  the  local  y  axis 

and  the  direction  in  which  the  force  is 
uesired 
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NOTATION  (Continued) 


t  -  direction  cosine  between  the  local  z  axis 

and  the  direction  in  which  the  force  is 
desired 

L  -  characteristic  length  of  the  problem  under 

consideration 


m 

n 

N 

q 

q 

R 

S 

t 


-  mass  of  a  molecule 

-  number  density  of  molecules  in  space 

-  total  number  of  molecules  that  impinge 
on  unit  area  per  unit  time 

-  mass  velocity  of  the  gas 

-  unit  vector  in  the  direction  of  the  mass 
velocity  vector 

-  gas  constant  for  a  particular  gas 

-  molecular  speed  ratio 

-  time 


T  -  temperature 

u,  v,  w  -  components  of  the  total  velocity  of  a 

molecule 


n,  v,  w  -  components  of  the  mass  velocity,  q,  of 

the  gas 

U,  V,  W  -  components  of  the  random  thermal 

motion  of  a  molecule 


x,  y,  z 


-  local  axis  system 


X,  Y,  Z 


-  body  axis  system 


<x 

C* 


-  angle  of  attack  defined  as  the  angle 
between  the  mass  velocity  and  the  longi¬ 
tudinal  axis  of  the  vehicle 

-  thermal  accommodation  coefficient 
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NOTATION  (Continued) 


z3 
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6 

ar 
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SUBSCRIPTS: 


r 


( i 


-  quantity  defined  as/3«  £p,j  ,  or  a 
general  geometric  angle 

-  semi-vertex  angle  of  a  cone 

-  direction  cosine  between  the  mass 
velocity  vector  and  the  local  x  axis 

-  direction  cosine  between  the  mass 
velocity  vector  and  the  local  y  axis 

-  direction  cosine  between  the  mass 
velocity  vector  and  the  local  z  axis 

-  one  half  the  included  angle  of  a  spherical 
segment 

-  collision  frequency  of  one  type  i  molecule 
with  type  j  molecules 

-  mean  free  path  of  a  molecule 

-  mass  density  of  the  gas 


-  incident  molecules 

-  reflected  molecules 

-  wall  condition 

-  an  underlined  symbol  denotes  a  vector 
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SECTION  I 
INTRODUCTION 


Analysis  of  the  trajectory  and  attitude  control  of  near  earth  satellites 
requires  a  knowledge  of  the  forces  and  moments  caused  by  molecular 
impingement.  A  detailed  analysis  of  these  molecular  impingement  forces 
and  moments  which  act  on  the  satellite  must  be  carried  out  and  incorporated 
into  the  vehicle  design. 

The  determination  of  the  aerodynamic  force  acting  on  a  satellite  is  a 
free  molecule  flow  problem.  The  first  extension  of  aerodynamic  considera¬ 
tions  to  free  molecule  flow  was  make  by  Tsien  (Ref.  1}  in  which  he  discussed 
the  method  of  defining  the  various  flow  regimes,  continuum,  transition  and 
free  molecule,  by  means  of  the  Knudsen  number.  Expressions  for  the  lift 
and  drag  coefficients  of  a  surface  element  were  derived.  Ashley  (Ref.  2) 
applied  free  molecule  theory  to  the  determination  of  lift  and  drag  coefficients 
for  several  shapes  and  discussed  the  possibility  of  aerodynamic  flight  in  the 
free  molecule  regime.  Detailed  studies  by  Stalder  and  associates,  (Ref.  3) 
developed  basic  free  molecule  flow  theory  and  applied  it  to  determine  lift 
and  drag  coefficients  for  several  shapes.  However,  certain  specific  assump¬ 
tions  regarding  the  temperature  (and  hence  the  velocity  and  momenta)  of  the 
re-emitted  molecules  which  are  not  applicable  in  satellite  environments  tend 
to  limit  the  practical  usefulness  of  these  results.  In  this  report,  any  perti¬ 
nent  value  of  the  temperature  of  the  reflected  molecules  can  be  used. 

A  closed  form  solution  for  the  free  molecule  aerodynamic  coefficients 
can  be  obtained  only  for  a  flat  plate  and  the  cylinder.  For  this  reason, 
several  approximate  methods  were  developed  by  Gustafson  (Ref.  4)  and 
Schrello  (Ref,  5).  Unfortunately  both  of  these  solutions  break  down  at  small 
molecular  speed  ratios,  and  at  small  angles  of  attack  for  all  values  of 
molecular  speed  ratio.  An  exact  analysis  is  necessary  in  many  cases  and 
is  presented  in  this  report. 
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The  interaction  of  molecules  with  each  other  and  with  container  walls 
falls  in  the  realm  of  the  kinetic  theory  of  gases.  The  fundamentals  of  this 
subject  are  contained  in  several  readily  available  texts.  However,  to 
facilitate  the  reader  without  a  background  in  this  field,  a  development  of 
the  basic  concepts  is  contained  in  Appendix  A  of  this  report. 
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SECTION  II 

FUNDAMENTALS  OF  FREE  MOLECULE  FLOW 

A.  ASSUMPTIONS  ON  WHICH  FREE  MOLECULE  FLOW  THEORY  RESTS 

Depending  on  the  governing  physical  phenomena  the  subject  of  fluid 
mechanics  can  be  broken  into  three  general  categories.  The  first  is  the 
ordinary  regime  of  continuum  fluid  mechanics  in  which  intermolecular 
collisions  are  the  important  physical  phenomena.  The  second  is  called  the 
"transition  regime"  in  which  molecule-molecule  and  molecule-surface  inter¬ 
actions  are  of  equal  importance.  The  third  is  the  free  molecule  regime  in 
which  molecule-surface  interactions  are  the  governing  phenomena.  The 
regime  in  which  a  particular  problem  falls  is  determined  by  a  consideration 
of  a  non-dimensional  parameter  called  the  Knudsen  number,  Kn.  This  is  de¬ 
fined  as  the  ratio  of  an  appropriate  mean  free  path(  A)  of  the  molecules  to  an 
appropriate  characteristic  dimension  (L)  of  problem  under  consideration, 
i.  e.  .KuH^.  A  detailed  discussion  of  the  determination  of  X  and  L  is  given  in 
Section  V. 

Based  on  Knudsen  number,  the  flow  regimes  are  arbitrarily  divided  as 
follows.  If  the  Knudsen  number  is  less  than  0.  01,  the  flow  is  in  the  con¬ 
tinuum  regime;  if  the  Knudsen  number  is  between  0.  01  and  10,  the  flow  is  in 
the  transition  regime;  and  if  the  Knudsen  number  is  grea.ter  than  10,  the  flow 
is  in  the  free  molecule  regime.  However,  there  is  some  experimental  evi¬ 
dence  which  seems  to  indicate  that  free  molecule  flow  may  occur  at  Knudsen 
numbers  as  low  as  3  in  certain  cases. 

That  the  Knudsen  number  is  actually  an  indication  of  the  flow  regime  of 
a  particular  problem  can  be  seen  by  the  following  consideration.  The  char¬ 
acteristic  length,  L,  is  determined  by  the  problem  of  interest.  If  the  Knudsen 
number  is  small,  the  distance  a  molecule  travels  on  the  average  before  it 
suffers  a  collision  with  another  molecule  is  small  compared  tc  the  geometry 
of  the  problem.  Since  wh  ere/^is  the  density  of  the  fluid,  this  implies 
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that  the  fluid  is  dense  and  thus  may  be  treated  as  a  continuum.  If,  however, 
the  Knudsen  number  is  large,  this  means  that  the  distance  a  molecule  travels 
between  collisions  with  other  molecules  is  large  compared  to  the  geometry 
of  the  problem;  and  >ince^W^^  this  implies  that  the  density  is  small  and  the 
gas  is  rarefied.  Thus,  the  fluid  may  no  longer  be  considered  a  continuum 
but  must  be  treated  as  an  assembly  of  discrete  particles.  For  intermediate 
values  of  the  Knudsen  number,  the  fluid  can  not  be  considered  to  be  a  con¬ 
tinuum  or  a  rarefied  gas.  The  flow  in  this  transition  regime  will  exhibit 
some  of  the  characteristics  of  both  continuum  and  rarefied  gas  flow. 

The  theory  for  predicting  the  aerodynamic  forces  on  a  body  in  continuum 
and  free  molecule  flow  is  available.  However,  at  the  present  time  there  is  no 
adequate  theory  for  predicting  the  aerodynamic  forces  on  a  body  in  the  transi¬ 
tion  regime. 

In  a  very  rarefied  gas,  which  is  the  main  concern  of  this  report,  the  mean 
free  path  of  the  molecules  is  much  greater  than  a  characteristic  body  dimen- 
sion.  The  basic  assumptions  of  free  molecule  flow  are:  1.  that  incident 
molecule-surface  collisions  are  much  more  numerous  than  incident  molecule- 
incident  molecule  collisions  and  2.  that  incident  molecule-surface  collisions 
are  much  more  numerous  than  incident  molecule-reflected  molecule  collisions, 
be. ,  the  incident  flow  is  undisturbed  by  the  presence  of  the  body  and  the  equi¬ 
librium  velocity  distribution  of  the  incident  molecules  is  changed  only  by 
collision  with  the  body.  Since  we  have  assumed  that  the  incident  and  re¬ 
flected  molecules  do  not  affect  each  other,  the  effect  on  the  surface  of  the 
incident  and  reflected  molecules  can  be  computed  separately.  Thus,  the 
flow  phenomena  are  governed  by  the  molecule-surface  interactions. 

In  order  to  compute  the  force  on  a  body  in  a  free  molecule  flow,  the 
molecule-surface  interaction  must  be  specified.  Most  surfaces  are  rough 
as  far  as  a  molecule  is  concerned.  The  impinging  molecule  hits  the  surface, 
bounces  around  in  the  spaces  between  the  molecules  or  atoms  of  the  surface 
and  then  at  some  later  time  leaves  the  surface.  This  interaction  is  usually 
considered  in  two  parts,  the  momentum  transferred  to  the  surface  and  the 
energy  transferred  to  the  surface.  The  momentum  which  is  imparted  to  the 
surface  depends  on  the  type  of  reflection  that  occurs.  There  are  two  types 
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of  reflection  that  may  occur.  They  are:  specular  in  which  the  molecule  hits 
the  surface  and  is  then  reflected  like  a  billiard  ball,  i.  e.  ,  the  angle  of  inci¬ 
dence  equals  the  angle  of  reflection  (Fig.  II-  1);  and  diffuse  in  which  the 
molecule  hits  the  surface  and  is  then  re-emitted  in  a  random  direction,  i.  e.  , 
there  are  as  many  molecules  with  a  particular  velocity  leaving  the  surface  in 
one  direction  as  there  are  leaving  the  surface  in  the  opposite  direction  and 
the  angle  of  incidence  of  a  molecule  is  in  no  way  related  to  its  angle  of  re¬ 
flection  (Fig.  II-2).  In  both  of  these  reflection  processes,  the  velocity  with 


Sr 


7777777777 


SPECULAR  REFLECTION  DIFFUSE  REFLECTION 
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Figure  ll-l 


Figure  11-2 


which  a  molecule  leaves  the  surface  depends  on  the  apaount  .of  energy  which 
the  molecule  transfers  to  the  surface  before  it  is  re-emitted.  The  degree  of 
equilibrium  attained  between  the  molecule  and  the  surface  before  the  molecule 
is  re-emitted  is  measured  by  the  energy  accommodation  coefficient,  C<  , 
which  is  defined  as 


0< 


where  E^  is  the  energy  carried  to  unit  area  of  the  surface  by  the  incident 
molecules,  E  is  the  energy  carried  away  from  the  unit  area  by  the  reflected 
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molecules  and  is  the  energy  the  reflected  molecules  would  carry  away 

from  the  surface  if  they  were  re-emitted  at  the  temperature  of  the  surface, 

Tw>  If  E^  =  E^,  there  is  no  energy  exchange  between  the  incident  molecules 

and  the  surface  andC(  -  0.  If  the  incident  molecules  and  the  surface  reach 

thermal  equilibrium  before  the  molecules  are  re-emitted,  E  =  E  and  0<=  1. 
n  r  w 

It  is  implicity  assumed  in  the  definition  ofO^  that  all  the  energies  associated 
with  those  molecular  degrees  of  freedom  which  enter  into  an  energy  exchange 
with  the  surface  are  accommodated  to  the  same  degree.  Experimental 
evidence  seems  to  indicate  that  this  is  true  for  the  translational  and  rotational 
degrees  of  freedom,  while  the  vibrational  degrees  of  freedom  are  practically 
not  affected  by  a  surface  collision.  If  necessary,  it  is  possible  to  introduce 
separate  accommodation  coefficients  for  each  degree  of  freedom. 

The  momentum  transferred  to  the  surface  is  specified  by  momentum 
accommodation  coefficients.  Maxwell  introduced  the  momentum  accommoda¬ 
tion  coefficient,  g,  defined  as  the  fraction  of  molecules  that  reflect  diffusely, 
the  remainder,  (/- g),  being  reflected  specularly.  At  the  present  time,  how¬ 
ever,  the  trend  is  toward  the  use  of  two  momentum  accommodation  coeffi¬ 
cients,  one  for  tangential  momentum,  {J~,  and  one  for  normal  momentum,  (J'/ ■ 
They  are  defined  as 

rr= 

r(  -  ru  v£ 
cr= 

f}  - 

where  is  the  tangential  momentum  carried  to  unit  area  of  the  surface  by 
the  incident  molecules,  Vr  is  the  tangential  momentum  carried  away  from  the 
unit  area  by  the  specularly  reflected  molecules,  rw  is  the  tangential  momen¬ 
tum  which  would  be  carried  away  from  the  surface  by  the  diffusely  reflected 
molecules  if  they  were  in  thermal  equilibrium  with  the  surface.  By  the 
definition  of  diffuse  reflection,  rw  is  equal  to  zero.  P^  is  the  normal  mo¬ 
mentum  carried  to  unit  area  of  the  surface  by  the  incident  molecules,  is 
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the  normal  momentum  carried  away  by  the  reflected  molecules  and  is  the 
normal  momentum  which  would  be  carried  away  by  the  diffusely  reflected 
molecules  if  they  were  in  thermal  equilibrium  with  the  surface. 

For  completely  diffuse  reflection,'^.  =  -0and  (T  -  1  regardless  of  the 
degree  of  thermal  accommodation.  If  any  specular  reflection  occurs,  the 
value  of  (T  depends  on  the  degree  of  thermal  accommodation  through  . 

If  there  is  complete  specular  reflection  and  no  thermal  accommodation, 
andO“  =0.  If  completely  diffuse  reflection  and  complete  thermal 
accommodation  occur,  C"  =  1;  if  completely  specular  reflection  and  no  thermal 
accommodation  occur,  cr  =  o.  For  any  type  of  reflection  between  these  two 
limits,  cr  depends  on  the  degree  of  thermal  accommodation  that  occurs. 

This  is  different  from£J"  which  is  equal  to  one  for  completely  diffuse  reflection 
regardless  of  the  degree  of  thermal  accommodation  which  occurs.  depends 
on  the  degree  of  thermal  accommodation  which  occurs  only  when  there  is  some 
specular  reflection  taking  place. 

Available  experimental  evidence  seems  to  indicate  that  for  the  type  of 
surfaces  used  in  most  satellites,  the  assumption  of  completely  diffuse  re¬ 
flection  is  legitimate  (Ref.  6).  This  assumption  has  been  used  in  the  devel¬ 
opment  presented  in  this  report.  However,  no  assumption  has  been  made 
concerning  the  temperature  of  the  re-emitted  molecules,  Ty.  Thus,  the 
results  presented  here  are  only  subject  to  the  assumption  of  completely 
diffuse  reflection. 


B.  DERIVATION  OF  THE  EQUATION  FOR  THE  FORCE  ON  AN  ELEMENT 
OF  AREA  IN  A  FREE  MOLECULE  FLOW 

We  now  compute  the  force  on  an  arbitrarily  oriented  element  of  ai  ea  in 
a  free  molecule  flow.  Since  it  has  been  assumed  that  the  incident  and  re¬ 
flected  molecules  do  not  interact  with  each  other,  we  first  compute  the  force 
due  to  the  incident  molecules  and  then  the  force  due  to  the  reflected  molecules 
and  add  them  together.  This  gives  the  total  force  on  the  element  of  area. 

There  is  one  additional  assumption  which  we  make  at  this  point  concerning  the 
geometry  of  the  surface  under  consideration,  i.  e.  ,  it  is  assumed  that  the  ele¬ 
ment  of  area  is  located  on  a  convex  surface.  If  we  do  not  make  this  assumption, 
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we  must  take  account  of  the  force  on  the  element  of  area  due  to  the  molecules 
that  have  been  reflected  from  other  parts  of  the  body. 

In  order  to  compute  the  force  which  the  molecules  exert  on  a  surface, 
it  is  necessary  to  know  the  velocity  of  each  of  the  molecules  in  the  gas.  We 
assume  that  both  the  incident  and  reflected  molecules  are  in  equilibrium  with 
themselves.  Kinetic  theory  then  states  that  their  velocities  are  specified 
by  the  Maxwellian  velocity  distribution  function  f.  The  velocity  distribution 
function  f  is  the  probability  that  the  velocity  of  a  molecule  selected  at  random 
will  be  in  the  range  u  to  u  +  du,  v  to  v  +  dv  and  w  to  w  +  dw,  where  u,  v  and 
w  are  the  components  of  the  total  velocity  of  the  molecule.  Then,  the  number 
of  molecules  in  unit  volume  of  space  that  have  velocities  in  the  range  u  to  u 
+  du,  v  to  v  +  dv,  w  to  w  I  dw  is 

A  f( Uj  ifj  ar)  du  d<J~  c(ur 

where  n  is  the  number  of  molecules  per  unit  volume.  A  discussion  and  der¬ 
ivation  of  the  Maxwellian  velocity  distribution  function  is  given  in  Appendix  A. 

If  a  gas  is  at  rest  macroscopically,  the  average  velocity  of  its  molecules 
in  any  direction  is  zero.  If  a  gas  has  a  macroscopic  velocity,  the  average 
velocity  of  its  molecules  in  three  mutually  perpendicular  directions  will  be 
the  components  in  those  directions  of  the  macroscopic  velocity  of  the  gas. 

The  macroscopic  velocity  of  the  gas  is  usually  referred  to  as  the  mass  velocity. 
Thus,  in  a  gas  which  has  a  mass  velocity,  the  mass  velocity  can  be  subtracted 
from  the  total  velocity  of  each  molecule.  The  resulting  velocity  is  called  the 
random  thermal  motion  of  the  molecule,  that  is, 

a  -zz  =  u 

U'  -  £r  —  y 
or  -ur  -  VV 


where  u,  v,  w  are  the  components  of  the  mass  velocity  and  U,  V,  W  are  the 
components  of  the  random  thermal  motion  of  the  molecule. 
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Consider  an  element  of  a  surface  which  is  fixed  with  respect  to  a 
body  axes  coordinate  system  X  Y  Z,  Figure  II- 3 .  The  coordinate  system 
X  Y  Z  is  stationary  with  the  gas  molecules  moving  relative  to  it  with  a 
mass  velocity,  q. 


We  now  chose  a  local  coordinate  system  x,  y,  z  on  the  element  of  surface 
area  with  positive  y  axis  being  the  inward  directed  normal  to  the  surface, 
the  x  and  z  axes  being  tangent  to  the  surface  element,  u,  v,  w  are  the  com¬ 
ponents  of  the  total  velocity  of  the  molecules  in  the  local  x,  y,  z  directions 
respectively. 

Consider  the  force  exerted  on  the  element  of  area  by  the  incident  mol¬ 
ecules.  Force  is  defined  as  the  rate  of  change  of  momentum.  We  have  assumed 
that  when  the  molecules  impinge  on  the  surface,  they  lose  all  their  momentum 
relative  to  the  surface.  Thus,  the  rate  of  change  of  momentum  of  the  incident 
molecules  is  obtained  by  integrating  over  all  possible  velocities  the  rate  at 
which  the  molecules  in  a  particular  velocity  range  transport  momentum  to  the 
surface. 

Only  molecules  with  a  velocity  component  normal  to  the  surface,  i.  e 
in  the  positive  y  direction,  will  hit  the  element  of  area,  The  number  of 
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molecules  in  the  velocity  range  u  to  u  +  du,  v  to  v  +  dv,  w  to  w  +  dw  which 
strike  the  element  of  area  in  time  t  are  contained  in  a  cylinder  of  height  vt, 
Figure  II-4. 


Figure  11-4 


The  number  of  molecules  in  the  velocity  range  u  to  u  +  du,  v  to  v^t-dv, 
w  to  w  4-  dw  per  unit  volume  is 

nL  -f  du.  dirdur  (0 

where  Hi  is  the  number  density  of  incident  molecules.  Therefore,  the  number 
of  molecules  in  the  velocity  range  u  to  u  +  du,  v  to  v  +  dv,  w  to  w  +  dw  that 
strike  unit  area  of  the  surface  in  time  t  is 

tr t  nc  fju dtrd u/-  [ z) 

The  number  of  molecules  striking  the  element  of  area  per  unit  time  is  ob¬ 
tained  by  dividing  equation  (2)  by  t.  The  number  of  molecules  in  the  velocity 
range  u  to  u  +  du,  v  to  v  +  dv,  w  to  w  +  dw  that  strike  unit  area  of  the  surface 
per  unit  time  is 

Hi  U~  ~f~ c/trL  efir  dur  (S; 

Each  molecule  transports  momentum  mu  in  the  x  direction,  mv  in  the 
y  direction  and  mw  in  the  z  direction  to  the  element  of  area.  The  force  is 
given  by  the  product  of  the  number  of  molecules  in  a  given  velocity  range 
that  hit  the  surface  per  unit  time  multiplied  by  the  momentum  they  carry  to 
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the  surface.  Thus,  the  force  in  the  x  direction  on  the  element  of  area  due  to 
molecules  in  the  velocity  range  u  to  u  +  du,  v  to  v  +  dv,  w  to  w  4-  dw  is 

mu.  n t  urf  Ju.  du-dur  (4-) 

The  total  force  in  the  x  direction  on  the  element  of  area  due  to  molecules 
in  all  velocity  ranges  is 


+00  +co  +00 


rn  r\i 

-00  -CD 


Litrf du-Jv-dus- 


(?) 


m  and  n^  have  been  factored  out  of  the  integral  because  m,  the  mass  of  each 
molecule,  is  a  constant  over  all  the  molecules  and  nt-  ,  the  number  density 
of  molecules,  is  assumed  to  be  uniform  locally  in  space,  f  is  the  Maxwell 
velocity  distribution  function  given  in  Appendix  A  as  ^ 

_L _ )  ^  2*rL 

ZTrRT/  c 

where  R  is  the  gas  constant  for  the  particular  gas  of  interest  and  u,  v",  w 
are  the  components  of  the  mass  velocity  ^  in  the  x,  y,  z  directions  respec¬ 
tively.  The  limits  of  the  v  integration  are  0  to  +  QC  because  only  molecules 
with  a  velocity  component  in  the  positive  y  direction  will  hit  the  surface. 
Molecules  with  velocity  component  in  the  y  direction  between  0  and  -  00  are 
moving  away  from  the  element  of  area  and  do  not  strike  it. 

A  similar  procedure  yields  the  expressions  for  the  force  in  the  y  and 
z  directions.  The  total  force  in  the  y  direction  on  the  element  of  area  due 
to  molecules  in  all  velocity  ranges  is 


+  00 


mnL 


+00 

r 


+oo 

r 

i 

ir 


■f  du  du-dur 


J  J 

-CD  0  -00 
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The  total  force  in  the  z  direction  on  the  element  of  area  due  to  molecules 
in  all  velocity  ranges  is 


+  CO  +  00  -hCO 

mn;j  I  urir~f ciu-  ckrd 

J  J  J 

-CO  0  -oo 


Ur 


The  vector  sum  of  eqs.  (5),  (7)  and  (8)  gives  the  total  force  on  the 
element  of  area  due  to  molecules  in  all  velocity  ranges.  However,  we  are 
interested  in  the  component  of  each  of  these  forces  in  a  particular  direction, 
for  example  the  normal  or  axial  force.  Let  ki  be  the  direction  cosine  between 
the  direction  in  which  the  force  is  desired  and  the  local  x  axis;  let  $  be  the 
direction  cosine  between  the  direction  in  which  the  force  is  desired  and  the 
local  y  axis;  and  let  t  be  the  direction  cosine  between  the  direction  in  which 
the  force  is  desired  and  the  local  z  axis.  Then,  the  component  in  a  particu¬ 
lar  direction  of  the  total  force  on  an  element  of  area  due  to  the  incident 
molecules  is  given  as 


Next  we  compute  the  force  on  the  element  of  area  due  to  the  re-emitted 
molecules.  Since  we  have  assumed  that  the  molecules  are  reflected  diffusely, 
they  are  re-emitted  from  the  surface  randomly  with  a  Maxwellian  velocity 
distribution  corresponding  to  some  temperature  Tr>  The  reflected  molecules 
are  thought  of  as  issuing  from  a  fictitious  gas  which  is  at  rest  relative  to  the 
surface  at  a  temperature  T  .  The  number  of  molecules  which  leave  the 
surface  per  unit  area  per  unit  time  is  the  same  as  the  number  of  molecules 
of  this  fictitious  gas  that  hit  unit  area  of  the  surface  per  unit  time.  Thus, 
the  number  of  molecules  leaving  the  element  of  area  per  unit  time  in  the 
velocity  range  uto  u  +  du,  v  to  v  +  dv,  w  to  w  +  dw  is 

nr(r{-r)'f'  (/o) 
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where  (-v)  is  used  because,  in  order  to  leave  the  surface,  a  molecule  must 
have  a  velocity  component  in  the  negative  y  direction,  n^  is  the  number 
density  of  the  fictitious  gas  from  which  the  reflected  molecules  are  assumed 
to  issue.  Since  force  is  defined  as  the  rate  of  change  of  momentum,  the  force 
on  an  element  of  area  due  to  the  reflected  molecules  is  given  as 

Force  =  (momentum  reflected).  ,  -  (Momentum  reflected) 

'initial  'final 


Since  the  reflected  molecules  are  initially  at  rest  on  the  surface,  (Momentum 
reflected)in.tial  =  0.  Thus,  the  force  in  the  x  direction  on  the  element  of  area 
due  to  re-emitted  molecules  in  all  velocity  ranges  is 


-mnr 


u  v 


u-f  d  u  dv-d 


CO 


-oO  -  0 0-00 


The  force  in  the  y  direction  on  the  element  of  area  due  to  re-emitted  molecules 
in  all  velocity  ranges  is 


+00  O 


-00-00 


(li) 


The  force  in  the  z  direction  on  the  element  of  area  due  to  re-emitted  mol¬ 
ecules  in  all  velocity  ranges  is 


+oo  o  +  00 


-mnr 


r 


Ur  c r 


-fd  U  du~dur 


I  J  0  J  J 

u  -CO-CO  -00 

The  vector  sum  of  eqs.  (11),  (12)  and  (13)  gives  the  total  force  on  the 
element  of  area  due  to  the  re-emitted  molecules.  However,  we  are  inter¬ 
ested  in  the  components  of  these  forces  in  a  particular  direction.  Using  the 
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direction  cosines  k,  1  and  t,  the  component  in  the  direction  of  interest  of 
the  total  force  on  an  element  of  area  due  to  the  re-emitted  molecules  is 


df^spLecTSD  _ 

7a  ~  mn> 


+c?  ?  f00 

?r  /  (f^Ui-J>tr+i:ur)u~TJuJtrcltV-‘  (V4) 


-00-00  -  o° 


The  component  in  a  particular  direction  of  the  total  force  on  an  element 
of  area  due  to  both  incident  and  reflected  molecules  is  given  by  the  sum  of 
eqs.  (9)  and  (14).  Thus 

+  CC+CO+CO 

=  /77nij  I  J (ku+JPtr  +  tuj-)(rf duclis-ciur 


-00  o  -00 


•hCO  0  -hCO 

r  ( 


-h  mnh 


(feu  +Jir  -htur)  li~  fdu du-dur 


05) 


-00  -0D  -00 


These  integrals  are  evaluated  in  Appendix  B.  It  should  be  noted  that  u,  v, 
and  w  are  the  total  velocity  components  of  the  molecules  while  the  distribu¬ 
tion  function,  f,  is  a  function  of  the  random  thermal  motion  of  the  molecules 
only.  Thus,  in  order  to  perform  the  integrations,  the  velocities  u,  v  and  w 
must  be  written  in  terms  of  the  mass  velocity  and  the  random  thermal  velocity 
of  the  molecules. 


Using  the  results  of  Appendix  B,  Eq.  (15)  becomes 


£Fj  _ 
dA  ~ 


( /  +  J Vi/ ~h  tn 


4  ft? 


{!€>) 
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where  y3  —  2.RT'  )  ~  2  RT  )  is  the  number 

density  of  the  fictitious  gas  from  which  the  reflected  molecules  are  assumed 
to  issue,  m  is  the  number  density  of  incident  molecules,  T\  is  the  tempera¬ 
ture  of  the  incident  molecules,  and  is  the  temperature  of  the  reflected 
molecules,  u,  ~,  w  are  the  components  of  the  mass  velocity,  q,  in  the 
local  x,  y,  z  directions  respectively. 

Equation  (16)  was  obtained  by  applying  Newton's  second  law  of  motion 
to  the  molecules  impinging  on  and  reflecting  from  the  surface.  We  now  need 
an  expression  for  nr,  the  fictitious  number  density  of  reflected  molecules. 

This  is  obtained  from  continuity  considerations.  Since  we  assume  that  mol¬ 
ecules  are  not  being  created  or  destroyed,  the  total  number  of  molecules  that 
impinge  on  an  element  of  area  per  unit  time  must  be  equal  to  the  total  number 
of  molecules  re-err.itted  from  the  surface  per  unit  time.  We  now  compute  these 
two  quantities,  set  them  equal  to  each  other,  and  solve  for  nf.  The  mechanics 
of  the  computation  are  carried  out  in  Appendix  B,  the  result  of  which  is 


nh 


[irJPtl  +  JvtfirJF )  +  jp 


(n) 


This  is  the  equation  for  the  component  in  a  particular  direction  of  the  total 
force  on  an  clement  of  area  in  a  free  molecule  flow  assuming  diffuse 
reflection. 
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At  this  point  it  is  convenient  to  introduce  a  dimensionless  quantity, 

S,  called  the  molecular  speed  ratio.  This  is  defined  as  the  ratio  of  the 

magnitude  of  the  mass  velocity  of  the  gas  to  the  most  probable  random  speed 

of  the  molecules,  c  ,  that  is 
m 


or  using  the  notation 


s  = 


where 


/ 

ZRTl  > 


C/n  =lfs  R'f? 


The  molecular  speed  ratio,  S,  is  a  measure  of  the  relative  importance  of 
the  mass  velocity  and  the  random  thermal  motion  of  the  molecules.  If  S  is 
small,  the  random  thermal  motion  is  important  in  determining  the  force  on 
an  element  of  area,  while  if  S  is  large,  the  effect  of  the  random  thermal 
motion  is  small. 

Let£_j^Y^  be  the  direction  cosines  between  the  local  x,  y,  z  axes 
and  the  mass  velocity  vector  cj_  respectively.  Then, 

“=  j  ^  =  *8  j  “'L  =  'T2 

fJF  =  r  iJF  =  y s 

Introducing  these  notations,  eq.  (18)  becomes  ^ 

fj  =  mni{w^k6  +J’r+t^[rf^’(l-h^r5)+lwe'!rs  1 

('+****)  +fa[*tfir(,+*0rs)  +  ^e~rs]J  ( (13) 
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Multiply  through  by  •  Then 


JE=  /? 

dA  2/3 


~JO- 

+j?jitfts)+se 

„  _ 


+4-(,+*tfn)+TM 


s)+e 


-rzsz 


(to) 


where  /?  ~tn  n: ■ 


Since  we  are  interested  ip  nondimensional  coefficients,  we  will  divide 
through  by  *£«*v 

Eq.  (20)  now  becomes 

dA  =  tv  f( S~”S] 

‘  '  ^  ^  2^.2 

jj  (2'^ 


2  2 

when  we  note  that  q  y3  =  S  .  This  is  the  equation  for  the  component  in 

a  particular  direction  of  the  total  force  on  an  element  of  area  assuming 

diffuse  reflection.  This  equation  is  exact  within  the  physical  assumptions 

of  kinetic  theory,  free  molecule  flow  and  diffuse  reflection. 

All  the  quantities  in  Eq.  (21)  are  known  except  for  T^,  the  temperature 

of  the  reflected  molecules.  This  can  be  determined  from  an  energy  balance 

for  the  surface.  To  simplify  the  computations,  T  will  not  be  specified. 

However,  T  will  be  assumed  to  be  a  constant  over  the  surface  under 
r 

consideration. 


ft 

W 
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To  remove  the  restriction  of  diffuse  reflection,  the  accommodation 
coefficients  as  defined  in  Section  II  A  are  introduced. 

cr  =  (tz) 

-ti 

cr'-  (23) 

*  -  £ 

Solving  these  expression  for"^  and  Pf*  gives 

Vr  -■  (l-tr)  Ti 
Pr=  (i-<r)p-  +  cr'Pw 

The  total  shear  force  on  an  element  of  area  is 

r  =  VL-zr  =  -0-<r)  r-  =  crrt 

The  total  normal  force  on  an  element  of  area  is 

P=  Pt+Ph  =  Pc  +f/-  <r')Pi  +  <r'%/ 

P  =  (Z-cr')P;  +  cr>w  (Zl) 


(2  4) 

(ZS) 

(Z£) 


Using  these  results  equation  (21)  becomes  j  j 

~dA  =  lU[rl6k  +  V)+(*-r’W]ftW”)+fc  e~rS] 

'  3  2  1 

+  1| Z^M.(i+^s)+  e  *  Jj  (ZB) 

-p 

An  alternate  form  of  the  expression  for  the  force  on  an  element  of  area 
for  the  case  of  diffuse  reflection  sometimes  appears  in  the  literature.  For 
diffuse  reflection  Cn=2jT  ™  unspecified,  and  EQ.  (28)  becomes 
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dA  =  ti  ^ /+^3j  +  # 

+  *{r£')J (1+****)+  f '£j2*pe/f.(i+*tfrs)  +  £  e/s|  (23j 


Comparing  equations  (21)  and  (29)  we  see  that  they  appear  different  even 
though  they  are  both  for  the  case  of  diffuse  reflection.  This  is  because, 
in  EQ.  (21),  use  is  made  of  the  fact  that  for  diffuse  reflection,  P  may 
be  computed  directly  instead  of  using  EQ.  (25)  to  obtain  P  .  EQ.  (21)  depends 
on  the  degree  of  thermal  accommodation  that  occurs  through  T  ,  the  tempera¬ 
ture  of  the  relected  molecules.  EQ.  (29)  depends  on£j"  which  in  turn  depends 
on  the  degree  of  thermal  accommodation  occurring. 
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SECTION  III 

DETERMINATION  OF  AERODYNAMIC  COEFFICIENTS  FOR  SIMPLE 
SHAPES  IN  A  FREE  MOLECULE  FLOW 

The  results  of  the  preceding  section  will  be  used  to  determine  the  nor¬ 
mal  force,  axial  force  and  moment  coefficients  for  several  different  bodies 
in  a  free  molecule  flow.  Since  the  variation  of  T^  is  a  complicated  and,  in 
general,  an  unknown  function  of  location  on  the  surface,  we  will  assume  that 
it  is  a  constant  over  the  surface. 

Since  a  body  in  a  free  molecule  flow  does  not  disturb  the  flow,  it  is 
permissible  to  resolve  a  complicated  shape  into  simple  parts  (cones,  cylin¬ 
ders,  etc.  )  and  to  compute  the  contributions  of  each  part  and  then  add  them 
together  to  obtain  the  total  coefficient  for  the  entire  vehicle.  Thus,  there  is 
no  need  to  consider  a  complicated  shape  in  general.  We  have  only  considered 
simple  shapes,  i.  e.  ,  flat  plates,  cylinder,  cone  and  spherical  segment,  be¬ 
cause  these  are  the  most  common  configurations  of  which  a  vehicle  is  com¬ 
posed.  One  must  be  careful  to  consider  the  shielding  of  one  part  by  another 
when  the  body  is  divided  into  simple  shapes.  The  equations  presented  here 
take  into  account  only  shielding  such  as  the  shielding  of  the  back  surface  of 
a  cylinder  by  its  front  surface,  and  thus  care  must  be  exercised  in  their  use. 
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A.  FLAT  PLATE,  BOTH  SIDES  EXPOSED  TO  THE  FLOW. 


Consider  a  flat  plate  oriented  as  follows. 


We  must  always  use  the  inward  normal  to  the  surface  under  consideration 
since  eq.  (21)  Section  II,  Part  B  was  derived  using  the  inward  normal.  From 
Section  II,  Part  B, 


-ras2 


dc  =  fyi  +tr/+>i-t)[rO+^)+  ^=s- e  J 


-h 


0+H*s) + V-+^^s)+  i*  e  *  S]J 


dA_ 

V 


0) 


The  unit  vector  in  the  direction  of  the  mass  velocity,  q,  is 


+  J  -^o( 
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Unit  vectors  in  the  directions  of  the  local  axis  systems  are  as  follows. 


For  the  bottom  surface 


For  top  surface 


■  x  =  cos « 
q  •  _y  =  sin  0< 
q  ■  %  -  0 


y  =  -  J 

z  =  -  k 

6  =  cosCtf 
y  -  -sinC< 
Yj  =  0 


Normal  F orce: 


For  the  normal  force,  the  direction  cosines  k,  1  and  t  are 

For  bottom  surface  For  top  surface 

k  =  x-j  =  0  k  =  0 

H  =  y  •  j  =  1  S.  =  -1 

t  =  z  ■  j  =  0  t  =  0 

Substituting  these  values  of  £j  .  tdj  t'  into  the  expression  for  dc  gives 

dC"  =  tt  e~*Vl+  &  ( 

0  fipr^r^m  .  .  \  l  .  •  .T  v /Tj  *n/Lrt^\ 


d° h  -  +  +  ~  «'  ]+ 1  (  W*) 

+■  fM[^(l+^rs)  +  75 

+  sjv^]  -& [>+"**$-  jjfpf 

where  the  relations  ^  =  -  %hottom  and  /tQp  =  -/bottQm  were  used 

to  express  everything  in  terms  of  (\  ->  . .  and  y.  At 

r  v  bottom  bottom 
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Tpj 

Assume/ I  -r 
becomes  jjOTTovv 


rL 

Tj 


Using  the  fact  that  erf  (-  x)»-erf  x, 


' rop 


(■ l£  I 

V  £  S  J 

Integrating  over  the  area  of  the  plate  gives 

CN 

i  CK 


2a*»,« 


?.  •  2 

—  S  4*. -r)  0{ 

e 


(z) 


Axial  F orce: 

For  the  axial  force,  the  direction  cosines  k,  %  and  t  are 
Bottom  surface  Top  surface 

k  =  x  •  i  =  1  k=l 

Z  =  y  .  i  =  0  J?  =  0 

l=2>«_i  =  0  t=0 

Expressing  everything  in  terms  °f£bottom>  ^bottom'  k  bottom'  dCA  18 

,2.  2 


^  "5  «* 

±ca  =^\f<*><x[r(i-t*ifrs)+  ^p.e  ]+'«»xfr(7+*y&tej)+^  JLxJ* 


Integrating  over  the  area  of  the  plate  gives 

CA  —  -4.  \z  c^tXA^CKAn/S^iCX  +  c  f 

A  4>jl  ^  J 


ZC^/X  ? 


(*) 
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PITCHING  MOMENT: 


The  pitching  moment  is  the  moment  about  the^I  axis.  Nose  up  moments 
are  positive.  Thus, 

dc"  =  t Hdc '* +  ft+x)Jciiv 

Substitute  in  for  dC^  and  dC/y  and  integrate.  This  gives 

C,,  =  P-CA  +  (L+  r) 


‘M 


Las/ 


Let 


-sltAJcx 


X^(SjOi)  =  2 


(S,o0  -  M^O( 


M 


(?) 


J 


Then  the  expressions  for  C^, 


C.  and  C..  become 
A  M 


CA 

CN 

CM 


J 
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If  the  edge  of  the  plate  to  which  L  is  measured  is  located  to  the  right 
of  the^  axis,  I  is  negative;  if  the  plate  is  located  below  the  2"  axis, 
hi  is  negative. 

B.  FLAT  PLATE,  ONE  SIDE  EXPOSED  TO  THE  FLOW 

We  will  now  consider  the  case  of  a  flat  plate  on  the  end  of  a  cylinder 
or  a  cone.  In  this  case,  one  surface  of  the  plate  is  inside  the  body  and 
thus  is  never  exposed  to  the  flow.  The  plate  is  located  as  follows. 


I 


Figure  111-2 

The  unit  vector  in  the  direction  of  the  mass  velocity  is 
q  =  i  cos  pi  +  j  si  n  of  ■ 

The  unit  vectors  in  the  direction  of  the  local  axis  system  are 

y  =  i*2.=  5s 

The  direction  cosines  £j  ^  are 

6  =  q_'X  =  -  sin£)(  ^  ^  '  J.  =  cos°^  3  *  -  =  ® 

Normal  F orce; 

For  the  normal  force,  the  direction  cosines  k,  1  and  t  are 

k* 1 = 

&  =  1  ‘  J  =  ° 

-t  -  2  •  j  ~  o 
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Substituting  into  the  expression  for  dc  gives 


A^cx^C(K(l+je^S^O(,)  +  e 


-5XQxLX 


]}« 


Integrating  over  the  area  of  the  plate  we  obtain 

a  C  Z  *  . 


St nr 


0) 


AXIAL  FORCE: 


For  the  axial  force,  the  direction  cosines  k,  1  and  t  are 

k  =  X-L  = ° 
l  =  1’L  =  < 


-f-  zr  *  L_  O 


Substituting  into  the  expression  for  dc  gives 

d CA  =  ^  £*««[<<*«  (/ + ^  e  s  w  j 

I  -jr  '  r-  /  7 

4-  (l-hjZAfSQr4  0<)  +  j~-C^c((l+^SC^oCj~l- 


Z  2  , 
-S^>o( 
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»  . 


1 

t 


Integrating  over  the  surface  we  obtain 


JV  [«T  «" ‘  ^  *]J 


Pitching  Moment: 


The  pitching  moment  is  the  moment  about  Lhe  axis.  Since  the  plate 
is  symmetric  about  the  “Z  axis,  the  axial  force  does  not  contribute  to  the 
pitching  moment.  Nose  up  moments  are  taken  as  positive.  Thus 


Integrating  wc  obtain 


C/YI  “ 


Let 


jz^CSjOl)  =  e 


?  Z 

—  see*  <X 


x'(Sjoi)  =  +  ^  +  c 


-  S  Vc< 


_  /SF 


2  5 


(I+jza£S<*4<x)  +  gjg  € 


7 
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Then  cN, 


C  .  and  C, ,  are 
A  M 


(10) 


J 


If  the  plate  is  to  the  right  of  the  J_  axis,  L  is  negative. 

C.  CIRCULAR  CYLINDER 

Next  we  will  consider  a  right  circular  cylinder  oriented  as  shown  below. 


*• -  YZ  - *■ 

* -  l  - * 

\ 

_ L 

J 

_ W 

\ 

\ 

\ 

« - Y  — 

l 

~  f  — ■ 

yr  \ 

/ 

w _ _ 

*■%  T 

x! 

V 

1  ^ 

(  z 

'C  BASE  OF  CYLINDI 

J  \ 

•R 

Figure  111-3 


We  will  neglect:  end  effects  in  computing  the  forces  on  the  cylinder.  Since 
the  cylinder  is  axially  symmetric,  it  will  always  be  possible  to  orient  the  body 
axis  system  so  that  the  velocity  vector  is  in  the  X-T  plane. 
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The  element  of  area  on  the  surface  of  the  cylinder  is 

d  A  =  r  djdy 

The  unit  vectors  in  the  direction  of  the  local  axis  system  are 

2S>  L 

V  =  jj  t&O  p  —  fe,  -dAn  p 
i  p  A  C&o  p 


The  unit  vector  in  the  direction  of  the  mass  velocity  is 

f  ~  At*>c<  +  _J  AAn  CX 

The  direction  cosines  6  ,  y  ,  are 

6  -  f  *  K  — 
ft  ~  f  • 

^  =  $  *  i  -  p 

The  integration  limits  for  ^  are  0  to  2^,  and  for  y  they  are  O  to  y^. 


Normal  force ; 

For  the  normal  force,  the  direction  cosines  k,  1  and  t  are 

k  -  X  •  j  =0 

r 

/  s  S  ’  J  = 

i 

~t  —  — '  -  =  4* 
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Substituting  these  expressions  for  JfJ  Yfj  l&j  Jtj  into  the  equation  for  dc 


dC/V  =  C  ^  -4*-n 

+  <£dL(i-WAfirs)  +  &iJ£  ji 


+  s  fir  e 


-  ■(  V- 


■s)+ite 


1 


Integrating  over  the  surface  of  the  cylinder  we  get 

2  7T  ZTT 

C^j  j  cl^  +  J J  -f  (f  d<f 

^  zsz{\  +  f  |T  [iJT^0<(l 

2  'jr'  \  *  z.tt  ^2.^1  -7 

■h  j  j  ~h  £Jzj  JfJ'i 

Simplifying  and  combining  terms  gives 

Cm  ~  j^(^nzo(  +  +  e  r  s d<i> 

+  /51 

JZ  as 


f  5  / 


The  integrals  are  evaluated  as  follows 

.2 V  .  ..  _  T*  .  ...  > 


C  CoA^SA^tsJ^  ~  Cirtj^fsdj1  -  4j  40^X5 
°  00 
re-^4^  fVv*-  ,/V^,ar^Xo(^  OB) 

X  X  Z 
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where  I  and  I,  are  modified  bessel  functions  of  the  first  kind  of  degree  zero 
o  1  b 

and  one  respectively.  Using  these  results  becomes 

+-/JE  (Zl_Ao^O<.[  (13) 

J  Z-  2  s 


AXLA  I,  FORCE: 

For  the  axial  force,  the  direction  cosines  k,  1  and  t  are 

k~  £‘  L  =  / 

i  -  l  *  i  =  o 

t  =  *  =  O 


Using  these  values  for  k,  1,  t, 

~  +sfrr  e  J yr-d^d'j 


Integrating  over  the  surface  we  obtain 


,  r2.Tf 

,  c^0<  r  ~2T^ 

+  — I  5  c 


5^ 
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Substituting  in  the  results  of  EQS.  (11)  and  (12)  gives 


j^2  SJ7F CertCX.  £  2  Cj  -f-  X 

-h  ZJEL  a.  e  ^  Zo(SAtr—)r 


Z  .  7  \ 


+■  (soo  o<  e 
S 


f/4) 


PITCHING  MOMENT: 


The  pitching  moment  is  the  moment  about  the  2  axis.  Moments  about 
the  ^  and  axes  are  zero  due  to  symmetry.  Nose  up  moments  are  taken 
as  positive.  Thus, 


dCM  = 


Integrating  over  the  surface  we  get 


C  =■  — 

A4  A.., 


Z  *  2 

+  Z,(*4z*)J  +  iff  e~^~X0  ) 

If''  ,  fz  ZlT  ZTT 

+J  rc  ffs  ^<x\~rzyz  wocj^  o{( £ 


LOCKHEED  AIRCRAFT  CORPORATION 


MISSILES  and  SPACE  DIVISION 


LMSC-448514 


Performing  the  indicated  integrations  gives 

f s )  e 
Aju4  iMif  gWg( 

+  A*,# e  s  Jl0( s- 


£  tin  <-tX 


'■) 


H  TfA^yy  0{  ^<>0{ 


(’*) 


Let 


c  Y  S  Aly\  ^  . 

X,  (s, a)  =  sJTm*,o<(z  ^  +  x\  e  «  [x,(£%*)  +X, 


5  2  Zg(. 

cY,  .  /_)v2  . 

I?(5;0<)=  HL^1<X 

d  2S 


9 


xf(h«)  =  ZslF^^«e  2  [x0(L^)+I,(^sf)j 


*W* 


_  err7  -  £  i  ■  t.  v 

-h  IE  Co^oi  £  2  (~£  °) 


Then  CLT,  C.  and  C,,  become 
N?  A  M 

,  ra.  r-r^Jn 


CN  ~  A 


{ x'  +/f 

\j 


X 


cr 


■ 


V 
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CM  = 


-  jlK.  rcY 
VW  3 


j 


If  the  point  about  which  the  moment  is  taken  is  to  the  left  of  the  base  of 
the  cylinder,  L  is  negative.  L  as  shown  in  the  figure  is  positive. 

D.  CONE  FRUSTUM 

We  will  now  compute  the  forces  and  moments  on  a  cone  frustum  of  circu¬ 
lar  cross  section  in  a  free  molecule  flow.  The  cone  frustum  is  oriented  as 
shown  below.  We  will  neglect  end  effects  on  the  cone. 


From  the  above  figure  wc  see.  that 

r  ~  y  iSin  & 
dr  —  i  dy 
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The  element  of  surface  area  is 

dA  -  rd$ds 

(ds)  —  (d'jj  +(cLr)  =  (d'))  =  (dy)  aj^ S 

ds  -  AZC-d cL'j 

ciA—  cLydity  ~  d^d<fi 

Cb-oS 


From  the  preceding  figure,  we  see  that  the  unit  vectors  in  the  direction  of  the 
local  axis  system  are 

X  =  -  j  <j>  -  h.  £&4<f) 

y  -  i  -f  j  —  _k  <j> 

2  —  L_  6  —  j  aA>i  &  C*C  (j>  -f  <j) 


The  unit  vector  in  the  direction  of  the  mass  velocity  is 

^  C*rC(X  -+*jJ  Of 


Since  the  cone  is  axially  symmetric,  it  will  always  be  possible  to  orient 
the  body  axis  system  so  that  the  velocity  vector  will  be  in  thejT-lTpl&ne  •  From 
the  definition  of  the  direction  cosines  6  ,  'If  ,  ,  we  find  that 

£  -  _X  ‘  %  =  —  o(  <j> 

—  C^oO,  CfO  £  —  Au^O^aI^A  C*b(f> 

» 
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Normal  Force: 

For  the  normal  force,  the  direction  cosines  k,  1  and  t  are 

k.  —  Jf  *  j  —  —  A^r\ 

J  =  V'i-  <+*»$<&>  4 


Substituting  into  eq.  (1)  for  k,  1,  t,  S  ,  gives  for  the  normal  force 

coefficient 


dcN  =  ir  J  -f  ^C^Oi (j)  0(  C&0~ ^  $ 

—  aJsy^O( <& Co-o  (l ^f5)  +  ske  ] 

(h*Y*s)  ,  ^T[jf(h^s)  *  &*)]ctA 


4.  £»-■» 

ZSz 


Simplifying  this  expression  we  get 

dCN = ^  *[*0+^) s] 

The  integration  limits  on  <f)  arc  0  lu  and  on  y  they  are  to  y Per¬ 

forming  the  indicated  integrations  wc  obtain 
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Investigation  of  the  above  integrals  gives  two  results.  First,  the  integrals 
involving  erf#" S  and  e  cannot  be  evaluated  in  closed  form,  and  second, 

these  integrals  are  symmetric  about  ^  =  Tf ,  i.e.  the  integral  from  0  to  I'TT 
is  equal  to  two  times  the  integral  from  0  to  Tf  •  Then  performing  the  integra¬ 


tions  that  can  be  accomplished  gives 

c  =  ^<4  ,  .  ,(*  > 

1T/Ty  j  2  d  +  Z 

4- Quad /2  )  I  qydpj^XSdj)  4-  f  S  df 

y4>  SJ^F1  J 

r — 1  V2  0 

/Tf  if  1  Z  1 

+  ^c<crttj^j>^tsd<t>J-+  Jc^^e^Su^Jr  0*) 


AXIAL  FORCE: 

For  the  axial  force,  the  direction  cosines  k,  1  and  t  are 

k  -  X  •  i  -  O 

H  ~  I'i  =  ^ 

-fc  —  £  '  _t  qm  X 
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Substituting  into  eq.  (1)  for  k,  1,  t,  £■  ,  y  ,  gives 

dCA  =  &n6  ydi  Cf  .  x,  % 

^  CuO  b  A  /  l\C*i0('Ll-ir'  o+i^o(c<rO&4j^&<&$+Co4o(CeJ'6 

^  t.  z 

+.  +  x  e-»  *j  j  ^ 


Integrating  over  the  surface  we  get 

C/,= 

+  0(zr  +f%ftSdt) + ^/J1 [Jf(f  Zt+fw*4! 


zT? 


sir..yV 


+  ft  e  ^ 


Simplifying  further  we  obtain 

/« 


j  Y^o. [T^<s +  f)  °Y2^  + 

r  <fr-  ^ 


I - 1  Vz  °  ZJ~^r'1'l> 

J7^  [^T~  ^QtA^jfjinjfirsdf 


^0(CpodJ '- C^j’JZsijfysdtfi^  -j-  ■4^-icf  J  g 


4-  <*LCr\ 
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PITCHING  MOMENT: 


The  pitching  moment  is  the  moment  about  the  axis.  Nose  up  moments 
are  positive.  Due  to  symmetry,  the  moments  about  the  2"  ?nd  yr  axes  are 
zero.  The  moment  due  to  the  force  on  an  element  of  area  is 


“ 


(L-l) 


— -  dtC. a 


Substitute  in  for  demand  dC^.  This  gives 

dCM  =  yVj 

.,-rr  r  r  ,  ,  /  J  , 


'  ]+jnr40+**frs) 

^Jf[!f0+M/rs)+  f.e'VjJ.Uut 


Performing  the  indicated  integrations  gives 

c"  ■  s?W  c-¥)]Hf* 

( zir  zir^  x  t.r 

+  J  Xsla$$S  d<j>  ■+  sjjp  J  i  <s^  ■A/tf'S' c£<f> 

0  ^  ^  ^ 
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J 


r 

Cc^Oi 


r 


vtt 


.vrr 


ftCerO#  c/j)  -f  j  tfCfrtcjjLyffr 5d<j> 


j  /■  2 .1Y  >  x  (  /  ZjTr  /MT 

sffF  J  c*y0cfl€  j  +  f  J  I-  J  c&o  f>  ■Qsyf frSdfj 

'  /-  fZTr  _?r^ 

/  tfcacfidf  ■+--£—  /  )f Qx><f>  (TS<d.<f>  +  ~z  <&4<f>C  dp 

Ja  s  4?  Jq  / 


_f_  Ain^S  1%  1  /Jit1 

2?  JtA 


Simplifying  and  combining  terms  gives 

1 _  da^S 


C  = 


{[h¥h¥: 


L 

7T 


Z/TTaA^  ca-4)Q(  -iU-vicf 
7T 


~h 


2  Of  Qk?oc  ^  fjg/i^ Xsdk  ~h  (^Z^ydu  ~h  S2n£ftS>cL§ 

Yi 


-f-  2  /a^v-i  o( 

S 


JtF  <x  c^aZS  +  (^*0(  4^n '• tsdff' 

.T  1  ~vzsz  ~\1 

+  tXC44& ( ?T £d<j> J  -j-  /  c&cjs  c  d<j>  J 

J  y*_  y3  j  r  fir 

-Z^nS  (--j-~)\ ^  ^ ^  j/^  ^5^ 

o  SJ7T  }a 

*  .  ,  r—,  rr 

-<2^>o  of  C^>  O-^Oyy  S  ~t.  ^'Tr  AjL*  <P(c*yO  d  y<U^  dj C~4 ScLfy 


4-  e  <^K)cy 

^  /~/~i3/2 


+•  ^  «  C*A^  CeJ-fciAjrsdifj  ~h  j j  (zo) 
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We  will  now  define  ihe  following  functions.  Let 


Co  ,  (■ 11 

^ )  —  2  “Tl* D(C&C 0( 4Usy\  S  ~j~  2  -dts** i  Of  Cer°  Of  (£j  (fSd-fy 

-h  CcyoS(zA^y)  cK  ~f  j C&xj)  JLn$ tfS dcj  ~h  f  £  *  Sc[(jl 

Xg  (SjCKjS)  —  J^jr-  ^>n  Ci.Co^'i 

+  Q&Z(f)  ^>Hr  Sdfi 

,CO  f ,  TT' 

X3  (SjO(;iJ  =■  TTc^OiM^aC^sS  +  A^n  f fZ<*0OC  + j C#0<f)J2tfysd<j> 

-h  Zc^oixU^  o {c^S(rc^)zf>M/^scL<j>-h  l&K  f&di*  ItS 

/  ,3/a  1  SM 

X4  ( oc ;6)  -  ^L-  ^noLc^o  SaI^  cf  +  JP? ^  Sp ’ca*0(^tjZ>$£A/te 


,-vzs\ 


4 -AU*o(C4K>(fJ  c+tyjtf}rsdsj  +  ‘ jco+St^dty 

co  .  y  0  ,'ir 

Xf  TT^Z  to<>0(.  -h  &  -j-  d^c(  ^ZboZo( -h  pzjj 


-\-  £  AyOy^  ct  Cf/O  l  j  OyO 


+mys^ 


X6  (SjOt;6)  (p.  +1^ 


~h  am-,  oi  Ce«f  /  c&J 


faytsdij  +  dtji 


where  Y  —  Ccm?o< dir>  6  -f  AOn  0<  c^cScgA^ 
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Using  these  relations,  the  expressions  for  C.  ,  ,  and  C..  for  the  cone 

■A.  JN  JV1 

become 


If  the  point  about  which  the  moment  is  being  taken  is  to  the  left  of  the  vertex 
of  the  frustum  of  the  cone,  L  is  negative. 

E.  SPHERICAL  SEGMENT 

We  will  now  compute  the  forces  and  moments  on  a  spherical  segment  in 
a  free  molecule  flow.  Since  the  spherical  segment  is  axially  symmetric,  it 
is  always  possible  to  consider  the  mass  velocity  vector  in  a  plane  of  symmetry 
We  will  neglect  the  effect  of  molecules  impinging  on  the  "back"  side  of  the 
spherical  segment.  The  spherical  segment  is  oriented  as  shown  below. 
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The  element  of  area  on  a  sphere  is 

A  A  —  e  A  © 

From  the  above  figure  we  see  that  the  unit  vector  in  the  direction  of  the  mass 
velocity  is 

=  L  <X 

The  unit  vectors  in  the  direction  of  the  local  axis  system  are 

X  ~  &  -  J  Cgrd  e  CoA  +  k  C&6 

V_  —  C&4  ©  4-  j'  ©  Cf-dfj)  —  /%  ^ 

x  can  be  arbitarily  chosen  in  the  plane  containing^  and  the  3C  axis  such 
that  it  makes  a  right  angle  with  y.  The  orientation  of  z_  is  then  fixed.  To  find 
the  orientation  of  z,  wc  use  the  fact  that  z  is  normal  to  both  x.  and  y  and  that 
the  dot  product  of  z  with  itself  is  one .  Thus,  let  Z  —  ~h  b  J_  ~h  C  ki 

Then 

2  «  X  =  O  =  CLdir\  Q  —  hoyO©  A-  A 

.  'j  —  o  —  CL  C&d  &-j-b  AMn  ©  C&ap  —  C  Ai* i  ©Aum  (ft 

/  =  a.z+  f+c* 

This  gives  us  three  equations  in  the  three  unknowns  a,  b,  and  c  which  we  can 
now  solve.  The  solution  of  this  system  of  equations  gives  for  z, 

v 

T&  —  A  Alstn  A  ~h  Jt  Cff-Ap 
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From  the  definitions  of  the  direction  cosines  ©  ,  we  find  that 

£  =  f'X  — 

Y~  f‘  V  =  <^&Cc46+Ain&Ainec&t$ 

‘  ?  —  dAsn  (X  <jjsr>  rf) 


The  integration  limits  on 


£ 


are  0  to  2  7?''  and  on  ©  they  are  0  to 


NORMAL  FORCE: 


For  the  normal  force,  the  direction  cosines  k,  1  and  t  are 

k  ~  X  *  J  ~  -QyJ&C&Xp 

Jl  -  V  *i  =  O  c^ap 

t  ~  2-‘i  -  p 


Substituting  into  Eq.  (1)  the  above  values  of  k,  1,  t,  &  ,  ft  ,  gives 


dcN  = 


r^[l  ■■y.c^sc^^ -f.q^cKc^e 

^  .  r  1 


ff 


+  OCoJ-j  -hA^n  O^^tiyd'pjJ  Y  ( I  +  ~h 

-h  ^2? 

a  .s8 


+ 


tp+Mjrs) + 

j jl  dj>d& 


1 


Simplifying  we  get 

dcH  =  4±  A; 


-rV!- 


dy  efv  ^  #  e~  J+ 


JUs>^  &  j 

a  ~ 


l 

v 


f  /"-^  f /+ +  s?  e" *  sl]}  d$dl 
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Integrating  the  preceding  equation  we  obtain 


c  =  X 
A 


■&,,  i  ir 

I  y  > : 


'  -o  wo 

<0/  r  ZTT 


J  J  &c/<pd&  J  $4in6JZA^tfS<&<l)c(e 

tv  fd,rzir  i  / ,%  rz7r 

+  J  d$de\+  ~zf  J  J'U'^G  c&4<f>d<pde 


6,^i7r 


l  L  ^ Z&  rsd<j>de)+  jJF [jMI  f j  '  j 

°  °  27r  J  \J.S  \  a  Jo 

■  if  ■  )  fe>  rzrr  rzs& 

I  j  ?f4oy\  O  c&t?(j>j2/i^}rsc(<t>a(G  J-f  r~2  /  / -<u^72<9  <^o<^  ^  c/<fid& 

J0  /  s  °  4>  Jo  J 


Investigations  of  the  resulting  integrals  shows  two  things.  First,  the 
integrals  involving  erf  J"5  and  e-0  cannot  be  evaluated  in  closed  form,  and 
second,  the  4  part  of  these  double  integrals  is  symmetric  about  =  nr  , 
that  is,  the  integral  over  (j?  from  0  to  Zlf  is  equal  to  two  times  the  integral 
over  ^  from  0  to  ~ff  .  Performing  the  indicated  integrations  that  yield  closed 
form  solutions  and  using  the  symmetry  property  of  the  <j>  integrals  gives 

r  _  r2  f  ,  ,  .  ye,rr 

/I  G,-\-Z  A<sr>0(eo4cx\ 

"v  L  "x? 

/\r'f1r  •  reir'1T  -yV 

+  (2*»fc  +  gijl  l^deq^xp  a^fysc(<)>cl&  l  Ajdnee.0 

y  0  O  ^  J7P1  Jg  Jq 

c<  (—Lf-!1  -  c^e,  4- 1/)  j.  J?r7 (cerftx J  J '■<*>& dty,z$  trsd^d* 

fe,rT  \  re'f'?r ■  °  °  -rV  1? 

+  Ainn<x\  j  ^o^Ji^^dejA-  -7J  dfidejr  (z?>) 


Jo  J0 
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AXIAL  FORCE: 

For  the  axial  force,  the  direction  cosines  k,  1  and  t  are 

k -  X'A  =  aLt)  e 
J  -  c^>e 

•fc  s  *  •  L  —  O 


Substituting  these  v.  ’uti,  of  k,  .  t  into  eq.  (1)  gives 

~  ^7T 1  •4x'r1  &  1  ( 0(  Cti-4  G  G  Co^d  +<4ri<X,C<f4Z0 

'V  L  .^S\ 

+  ajl^VL/u^ &c^ec#4<i>)ly(l+aY)rs)-h  rj=?f  j 

+  ^(,+^s)+fz  e  *V]J^ 


Simplifying  the  above  expression  and  integrating  over  the  surface  gives 


£4  =  J<&< ><x£  J  J  fa i„  +  J  Jfa**oj>4firsdtyct6 

j_ f'6lfz T.  ~yzs'l  1  /’re,fZ7r 

+  J  l  dj>dsj  +  —ef  J  J ec&tedfyde 

reirZ7r  \ 

-h  I  AineC0rSQ&l{VScL$dQ] 


J 


-M*. 

m 


3-Z7 


LOCKHEED  AIRCRAFT  CORPORATION 


MISSILES  and  SPACE  DIVISION 


LMSC-4485 14 


Performing  the  indicated  integrations  gives 

S'  —  r _ ^  2  a  i  'X'-  - ; 


C-  rlC 

—  f  rr^ZC<^Oy^S>,+  Jt-  4X^0,  +(2caSLcX-h-L  )  tsd^de 

Am6  i  25  {  s  JJo  Jo  A  ~ 

+- 2  ^0(0^0^ J  J  fj^rj  J^ee 

~~ — c<^ac<(l~'  ^20/J  X$d<i>dG 

f6i  rTT  \  r6‘ r7r  ~&ZSZ  ~fl 

4 'jUmtxJ  j  ^ze^ec^^rsd<fid&J+£J  jM^ec^ee  d$de jjC^J 


3  S 
e,  rir 


PITCHING  MOMENT: 

The  pitching  moment  is  the  moment  about  the  Z  axiB,  Due  to  symmetry, 
the  moments  about  the  X  and  Y  axes  are  zero.  Nose  up  moments  are  positive. 
The  pitching  moment  is  given  by  the  following  equation. 

oLCm  -  J1—.  (l  trees')  die N  -  d-  r*s»ei#t>$<LcA 


Expanding  this  expression  we  obtain 

Cai= +vv  hii  ej}^^e^d9 


ert  .*•"  royt"  _  ^  jJ'* 

+  [  s  c&g  XsjLjdle  +  ^h=;  l  A^ectoee  d<$de 

'h  Jo  ~  17  v 

1*11  Q)  riff 

*-&(l  &L->£  SC<y*  9  asofyurtf , Kzdtf. 


I  j  C64&C4<$d$<l&+  j 


Arzir 


&,zr 


?[2S 


j  &  (-9^  9  Co4$d<f)cl&  "i~J  j  ^ac-A&Co4&C^s^ji/^ ^sd<pd&J 
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+  ^ 


e,  r  Z7r 


-b  Jo 


■Adecco e <&>j>  e  d$d@ 

rr^rzlf  . 


,6,^ 


r  r  r^1  r  r'r 

,  ,e,rz?r  _^sz  ~j  r6'fz7r 

+  S~fiF'  I  J  4^y'Z&c*o<t>s  ^ y  2S*  ( J  J'^Z0cM®c*<)<i,d<pd9 

r&,  rZIr  \ 

+  I  A<sr\  ©  C&4  &  C&d  <j}  yS  d(j> did  J 

4  4  / 

P _ j  g  2  'jr* 

^J-y  j^~- '(  j  ^  X<d^zec^ec*c<f>d<t>de 

r&ifz7r  re'rZ7r  <;  e  ~i? 

~l~  J  Ysd.^)(i& -f-  j— -ZJ  J ^dz&cgt<oec<^>^  e  ^  sd(j>dy 


Simplifying  the  above  equation  and  performing  the  integrations  that  yield 
closed  form  solutions  gives 

I—  d  C  ? 

CM  =  y~  CN+-~ - jrr^o<c^)c x(c**erY>Ji8l) 


L.  C  \ 

M—  j-~  £*+“ —  \ rr^ a c^oo(  fa e, - c&a 3e, ) 

H  C&QcJt  Adsr\  <x^  j  'J  Ai^ecod'e  jyyvSdL^de—^  j‘dyd&ci&Z<fijZrtrfXSdfe{< 

1 6  r  W 

4 •  z(a<^zc( -  Qr^cd)  (  CQ~d&C0-o§je/itf  XSd<t>d& 

•4  4>  z  z 
P  /  r6>rlr .  4  s . 

't  ~j_r1  ^  Y\  ay-'>i  &<<&&  s  d§d& 


r*  rv  .  z  ,  -*V 

—  <dXlO(  AOy\  ecthfep  e  dfdB 
Ja  J* 


(2  5-j 
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Now  let  us  introduce  the  following  notations, 


-X,  (SJ°()&')  ~  V*'4^C<CrtCX<UsylZ&/  +  2*^\0(Gx)C(f  ‘f^QCtrtQJ^tfsdjxie  ^ 

-  £  r  Jr  -*0  Ja 

+■(2.  Ainrfa(-h  ~ | Jj  lAOydec^^Jl^'Xsd'fide 

'  VI  ~V5 

'  diddle 


JC  JC5 

i _  ex*iv»c<  ('e,{rr.  -a'V'- 


^2?^^')  ®  &!-4*»x(Ss£2t-e*e,  +  £j 


+ML 


0\rTr 


(<*>«('( 
'  >0  4j 


&C&4 

Jlr+ii 


i  /- 6i  „  .y^e® 

+  ’si\  I  **”Zec*°$e  ct<t>de 
~o  4> 


X*3  (S^CXj  ©,)  =  ^2C^Z0(+ 

r6,rT  °  a  f 6‘fr 

-bS/u~uxc04o<l  l^z6c&<l)j24tfysdt<f>d&+!z£^\  b^ee  d<f>c(e 

•to  Jo  sJ^n[  J 

,  3/2  » _  firlf  *  0 

x\  ( Sj  <X,  ©,)  =  fg-  <*»*(/- anty+Kl  (wj  JaU> &  cj-e  w/ndbcte 

•  fe,rr-  »  \  W  ~ys\z 

J  ^  jj^vsc^ee  d<j>d& 

SP 

JZc  (Sj  <X .  ©J  =■  7T-Xv)  c<  ©,  -  <^o3©,  j 


'  0  *T  '  /  ^  ^ 

+  j^d^Bcoc’QM^yzd^d^  —  j  {/u#?6 ceJ'fywtd tfSdjide 'j 

0  6,e,rr  *  o 


o  o 


+  2  j  Lm^Z0  cc4eCe*}<})&ttfXSd<t!dt& 

Jo  J0  ,  , 

2.  /  re,rr  ~*V  r*'rr  -tszsz  \ 

€  ddd&-c^o(^J^Z9Q^<f>e  dfidej 


J 


where  —  <^’>d  Qyd  &  ~h^y i  (X  &  QK>(^ 
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If  the  point  about  which  the  moment  is  taken  is  to  the  left  of  the  center  of  the 
sphere,  L  is  negative . 

The  foregoing  analysis  resulted  in  closed  form  solutions  for  the  forces 
and  moments  on  a  flat  plate  and  a  cylinder.  However,  it  was  impossible  to 
obtain  analytic  expressions  for  the  cone  and  spherical  segment.  Since  there 
are  several  approximate  solutions  available  (4),  (5),  it  was  felt  that  another 
approximate  solution  would  not  serve  a  useful  purpose.  Further,  both  approx¬ 
imate  solutions  fail  at  low  molecular  speed  ratios,  and  at  small  angles  of 
attack  for  all  values  of  the  molecular  speed  ratio.  Thus,  it  was  decided  to 
evaluate  the  exact  expressions  numerically  so  that  the  approximate  solutions 
could  be  compared  with  the  exact  solution. 

All  of  the  expressions  for  the  force  and  moment  coefficients  of  the  shapes 
investigated  were  able  to  be  written  in  the  following  functional  form, 

C  =?(hn -)\i,  (w,jl 

where  gib1-)  is  a  function  of  the  length  dimensions  of  the  shape  only, 
and  (  Sj  dj  j3  )  are  functions  only  of  the  molecular  speed 
ratio,  J  ,  angle  of  attack,  ,  and  in  the  case  of  the  cone  and  spherical 
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segment  a  geometric  angle,  j3  .  For  the  cone,  the  geometric  angle  /3  is 
the  semi-vertex  angle,  ,  and  for  the  spherical  segment,  /3  is  half  the 
included  angle  of  the  segment,  ©j.  Since  the  I's  are  functions  of  S  ,  Of  , 
and  only,  they  may  be  considered  as  general  free  molecule  flow  functions 
which  are  independent  of  the  length,  radius,  etc.  of  the  particular  shape  under 
consideration.  Once  they  are  known,  the  aerodynamic  coefficients  for  the 
particular  shape  can  be  computed.  In  general,  the  I's  were  seen  to  be  com¬ 
plicated  integrals  which  can  be  evaluated  in  closed  form  only  for  the  flat  plate 
and  cylinder.  For  this  reason,  they  were  evaluated  on  the  IBM  7090  computer 
and  are  presented  in  tabular  form  for  the  following  ranges  of  the  parameters 
S  ,  CX  ,  6  ,  and O  ^ :  I£S±/3j  O^CXSlSo'j  S°£  0^/80° 

where  6  and  ©  ^  increase  in  increments  of  5°, 

The  0(  range  was  chosen  to  be  0  to  ^because  the  magnitude  of  the  I's 
in  the  range  of  Of  from  0  to  •  IT  is  the  same  as  in  the  range  0  to  TT •  It 
was  felt  that  values  of  5  from  1  to  13  would  cover  the  range  of  molecular 
speed  ratios  likely  to  be  encounted  in  practice.  The  use  of  the  tables  of 
the  I's  for  computing  the  free  molecule  aerodynamic  coefficients  of  a 
particular  body  will  be  discussed  in  Section  V. 


1 

I 
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SECTION  IV 

COMPARISON  OF  THE  EXACT  AND  APPROXIMATE 
FREE  MOLECULE  THEORIES 

Two  approximate  theories  (refs.  4  and  S)  for  determining  free  molecule 
aerodynamic  coefficients  have  been  presented  in  the  literature.  Gustafson's 
theory  (ref.  4)  consists  of  an  approximation  of  the  physical  phenomena  which 
occur  in  a  free  molecule  flow.  In  this  theory  it  is  assumed  that  the  random 
thermal  motion  of  the  molecules  can  be  neglected  in  comparison  with  their 
mass  motion.  This  physical  assumption  results  in  a  infinite  molecular  speed 
ratio,  S  =  00  and  is  known  as  the  hyperthermal  free  molecule  theory.  The 
kinetic  thecry  should  give  the  same  results  as  this  hyperthermal  theory  when 
S  is  large.  In  the  kinetic  theory  free  molecule  expressions,  S  is  always  mul¬ 
tiplied  by  a  function  of  the  angle  of  attack  which  goes  to  zero  astX"*1  0. 
Thus,  even  for  large  S,  there  will  be  a  discrepancy  between  this  approximate 
theory  and  the  kinetic  theory  approach  for  small  angles  of  attack. 

In  the  hyperthermal  theory,  the  incident  molecules,  which  have  mass 
motion  only,  impinge  on  the  surface  and  give  up  all  their  momentum  to  it. 
They  are  then  assumed  to  be  reemitted  diffusely  with  a  temperature  corre¬ 
sponding  to  the  wall  temperature  T  . 

For  a  convex  body,  the  force  on  an  clement  of  area  in  this  type  of  flow 
is  given  as  follows.  By  Newton's  second  law  of  motion,  the  force  is  equal  to 
the  rate  of  change  of  momentum. 


Figure  IV-1 
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x,  y,  z  is  a  local  coordinate  system  on  the  element  of  area  with  y  in  the 
direction  of  the  inward  normal  to  the  surface.  Then 


dF  =  pfXdA 


where  y  is  the  direction  cosine  between  the  mass  velocity^  q^nd  the  inward 
normal  to  the  element  of  area,  y.  Nondimensionalizing  we  obtain  for  the 
force  coefficient 


dc  =  zY 


£a_ 

A y 


(0 


Since  the  molecules  have  mass  motion  only,  the  total  force  due  to  the 
incident  molecules  is  in  the  direction  of  the  velocity  ^  which  by  definition 
is  the  drag.  Thus 


dc0  =  ay  d&-  (2) 

If  we  take  the  projection  of  A  normal  to  the  velocity  as  the  reference 


area,  then  becomes 


(3; 


based  on  the  projected  area  normal  to  the  velocity  q. 


Since  the  total  force  due  to  the  incident  molecules  is  drag,  we  see  that 
there  is  no  component  of  force  normal  to  the  drag,  i.e.  there  is  no  lift  due 
to  the  incident  molecules.  This  is  true  only  for  this  theory  which  is  based 
on  tie  assumption  of  no  random  thermal  motion  of  the  molecules. 

To  find  the  force  on  the  element  of  area  due  to  the  reflected  molecules, 
we  take  the  limit  as  S— >  OO  of  the  kinetic  theory  expression  for  the  force  on 
an  element  of  area  due  to  diffusely  reflected  molecules.  From  part  8  of 

IS* 

Vi  section  II, 
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^ -  ij-f  [*spr”(i+e«/rs)  +  er 5  ]dA  (4) 

Taking  the  limit  of  this  expression  as  S  — *■  CC  gives 

J_F  =  FF  SXpj™  cLA 

‘nerLecTev  g  p  J  j\ 

Nondimensionalizing  gives 

dCriePLecTB-e  _  JtP )f  f-fw  1 

This  force  is  normal  to  the  surface  because  the  diffusely  reflected  molecules 
have  a  net  charge  in  momentum  normal  to  the  surface  only.  The  shear  com¬ 
ponents  cancel  each  other  out.  To  find  the  total  force  in  a  particular  direc¬ 
tion  due  to  both  the  incident  and  reflected  molecules,  Eqs.  (2)  and  (6)  must 
each  be  multiplied  by  the  appropriate  direction  cosine  and  then  added  together . 
The  total  force  on  a  body  in  a  particular  direction  is  then  obtained  by  integrat¬ 
ing  the  total  force  equation  over  the  surface  of  the  body.  The  results  for  a 
cylinder,  cone  and  spherical  segment  are  presented  in  reference  (4). 

The  approximation  by  Schrello  (ref.  5)  is  a  mathematical  assumption  which 
was  introduced  in  the  kinetic  theory  expressions  for  the  force  and  moment 
coefficients  in  order  to  obtain  a  closed  form  result.  This  approach  makes 
no  assumptions  about  the  physical  processes  occurring  in  a  free  molecule 
flow.  The  mathematical  assumption  does  result  in  a  criteria  which  a  body 
must  fulfill  before  the  results  can  be  used  to  compute  the  force  and  moment 
coefficients.  This  criteria  is  that  the  average  normal  component  of  the  molec¬ 
ular  speed  ratio  over  the  exposed  body  surface  must  be  greater  than  or  equal 
to  unity.  This  approximation  also  fails  at  small  molecular  speed  ratios,  and 
at  small  angles  of  attack  fur  all  values  of  the  molecular  speed  ratio. 
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The  exact  theory  includes  co  .'.-ion  :>f 

whereas  the  approximate  tlieori  ■  •  * 

elusion  of  the  random  motions  \  u 

number  of  molecules  impinging 
This  is  illustrated  by  considering  the  au 
molecules  impinging  on  unit  area  per  unit  time, 
and  N.  as  given  by  the  theory  of  reference  4.  From  App  -ud'i  1 
theory  gives  IT,  when  the  molecular  speed  ratio  is  introduced, 


l>"  •  •  -!■■»• 

i  .  .  i..  iv.  An  v 
. .)  be  the  exact  thco  •• 
'  the  exact 


Ni=r'‘l[{(i+«JK)+  i%F 

According  to  the  hy perthermal  theory  (ref.  4), 

Ml  =  hi  fit 


-jr*sz 


(V 


Then. 


^  -  ~ =  I  +  zwf 


2  2 
-Jf  52 


(9) 


In  Figure]y-2j  is  plotted  against  the  angle  of  attack,  Of  >  of  the  surface 
for  several  values  of  the  molecular  speed  ratio,  S.  As  seen  in  Figure  JE-—  £  , 
for  values  of  S  greater  than  5,  which  are  pertinent  to  a  near  earth  satellite, 
the  random  motion  effect  is  confined  to  small  surface  inclinations.  As  such, 
the  primary  effect  would  be  an  increase  in  the  axial  force  coefficient,  with  a 
relatively  small  effect  on  the  normal  force  coefficient.  This  is  shown  in 
Figure  IV- 3  where  the  axial  and  normal  force  coefficients  for  a  surface  ele¬ 
ment  are  plotted  as  a  function  of  surface  inclination.  From  Figure  IV- 3  we 
can  readily  see  that  the  approximate  theory  is  quite  adequate  for  the  prediction 
of  the  normal  force  at  all  surface  inclination  angles,  but  is  inadequate  for  the 
prediction  ol  the  axial  force  when  the  surface  is  closely  aligned  to  the  free 
stream . 
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It  is  interesting  to  apply  these  results  to  an  idealized  satellite  shape  such 
as  that  shown  in  Figure  IV-4.  The  aerodynamic  force  and  moment  coefficients 
are  presented  in  Figure  IV- 5.  The  computations  were  performed  for  a  mo¬ 
lecular  speed  ratio  of  7  and,  assuming  complete  thermal  accommodation, 
fv*'  =  0.425  which  are  typical  of  the  conditions  at  satellite  altitudes, 
i.  computing  the  coefficients,  the  forward  end  of  the  cone  and  aft  end  of  the 
cy]  n  ••  -ore  considered  as  flat  plates.  The  results  using  the  theory  of  refer- 
er.o  5  ate  J  for  clarity  of  presentation  since  they  would  fall  between 

the  exact  anti  ro\;  :i«  .  results  of  reference  4. 

A">  se.'n  in  F i ■  i  ■  I',  •*,  th  •'’•iinary  difference  between  the  approximate 
ami  e  \  M  m..  '  •••"  Iws  :i.  •  he  prcchc<  u>;>  •f  the  axial  force  coefficient  at  small 
ts'ea  of  .  'iilferenc,  •  as  u-o  v*  previously  seen,  to  the 

•'  the  I'.--  "'om  ■  so .  t  ...  .ions  .  the  cun.t  ?l  surface  elements 


he  0'v£li.  t  -.L  I'*,  v.i  ij .  row 


■n  in  F i •  i 1 


1  v  :i  of  ;  i  i,i  . 


■  i  i  1 1  c  r  o  n  c 


■ •  h'l  r 


m-Wo  p red ■  i  t ion s 

face  of  the  l  \  ’  i  ndi 


•mi  'he  v  yle.ui  •  ace  elements 

<•  di!!' i  •  once  !•••  tv. -a  exact  and  approx - 
i  ;<•  mgement  on  the  sur- 

'  .s  of  the  results  shown 

-u.  is  a  linear  function 
<n  uu  ‘-o\  ,  ited  in  Figure  IV-6. 

..he  approximate  theory 
of  CN  ..  all  angles  of  attack,  but 


in  Figure  IV  2.  1,.  .•..•.•jiiiimie  o  •  n,  is  a  linear  function 

of  the  ^/O  of  the  cylindric_  <n  vnt  .  o\  ,  ited  in  Figure  IV-6. 

As  noted  from  the  analysis  of  a  urn.  rfu<  ..he  approximate  theory 

appears  adequate  for  the  prediction  of  ..  all  angles  of  attack,  but 

is  inadequate  for  the  prediction  of  for  an  clong<_  cd  body  at  small  angles 
of  attack. 

The  effect  of  surface  inclination  on  tlse  agreement  between  the  exact  and 
approximate  theories  is  illustrated  further  in  Figure  IV-7  in  which  at 
0(  =  0  is  plotted  against  L/d  for  a  15°  and  a  5°  cone,  and  a  cylinder. 

End  effects]  and  the  contribution  of  the  reflected  molecules  have  been  neglected 
so  that  the  discrepancy  between  the  curves  is  due  entirely  to  the  incident  mole¬ 
cules  • 

In  Figure  IV-8,  C^j  and  for  a  spherical  segment  are  presented  versus 
angle  of  attack.  From  the  preceding  discussion,  the  close  agreement  between 
the  exact  and  approximate  theories  for  this  configuration  is  to  be  expected. 
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Flgur e  IV- 2  Difference  In  the  Number  of  Molecules  Impinging  on  Unit  Area  Per  Unit  Time 
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0  30  60  90 

ANGLE  OF  ATTACKS 


Figure  IV-3  Aerodynamic  Coefficients  for  a  Flat  Plate  with  One  Side  Exposed 
to  the  Flow  Assuming  Complete  Thermal  Accommodation 


A"f  n  Plan  Farm  /.r»o  of  the  Plata 
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Figure  fV-5  Aerodynamic  Coefficients  for  Idealized  Satellite  Shape 
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Figure  IV -8  Aerodynamic  Coefficients  for  a  Spherical  Segment 


LMSC  448514 


SECTION  V 

USE  OF  THE  FREE  MOLECULE  FLOW  TABLES  FOR 
DETERMINING  THE  AERODYNAMIC  COEFFICIENTS 
OF  A  PARTICULAR  VEHICLE 


The  purpose  of  this  section  is  to  discuss  some  of  the  difficulties  which 
might  be  encountered  in  analyzing  a  vehicle  in  a  free  molecule  flow. 

A.  DETERMINATION  OF  THE  KNUDSEN  NUMBER 


The  first  step  in  analyzing  a  vehicle  is  to  determine  if  a  free  molecule 
flow  exists.  This  question  is  answered  by  computing  the  Knudsen  number, 

Vl  ,  where  \  is  an  appropriate  mean  free  path  of  the  molecules 
a. id  L  is  an  appropriate  dimension  of  the  vehicle.  If  Kn  is  greater  than 
10  (Kn>  10),  a  free  molecule  flow  exists.  For  purposes  of  dete rmining 
whether  or  not  an  entire  vehicle  is  in  a  free  molecule  flow,  the  largest  di¬ 
mension  of  the  vehicle  is  the  appropriate  L  to  use  in  computing  the  Knudsen 
number . 

The  determination  of  the  appropriate  mean  free  path  is  more  difficult. 

(A  discussion  similar  to  the  following  is  given  in  reference  7.)  The  magnitude 
of  the  mean  free  path  depends  on  the  coordinate  system  in  which  the  motion 
of  the  molecules  is  observed.  Since  wc  arc  interested  in  a  body  which  is  mov¬ 
ing  with  constant  velocity  through  a  gas  at  rest,  there  are  two  coordinate 
systems  with  respect  to  which  the  motion  may  be  observed.  We  may  refer 
the  motion  to  a  coordinate  system  that  is  stationary  with  respect  to  the  gas  or 
we  may  use  a  coordinate  system  that  is  fixed  in  the  body.  Since  the  expres¬ 
sions  lor  the  aerodynamic  coefficients  have  been  written  with  respect  to  a  co¬ 
ordinate  system  fixed  in  the  body,  this  is  the  coordinate  system  that  will  be 
used  to  Compute  the  mean  free  patli  of  the  molecules. 

There  are  four  mean  free  paths  that  can  be  defined:  The  mean  free  path 
of  the  incident  molecules  relative  to  themselves,  \ci.  >  the  mean  free  path 
of  the  incident  molecules  relative  lu  the  reflected  molecules,  >  the  mean 
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free  path  of  the  reflected  molecules  relative  to  the  incident  molecules,  Ari  , 
and  the  mean  free  path  of  the  reflected  molecules  relative  to  themselves , 

Kr  ’ 

Since  the  Knudsen  number  is  used  to  determine  when  a  free  molecule 
flow  exists,  it  must  indicate  when  the  assumptions  on  which  free  molecule 
flow  is  based  are  satisfied.  These  assumptions  are:  1.  it  is  assumed  that 
molecule  -  surface  collisions  are  the  governing  physical  phenomena,  i.e.  in¬ 
cident  molecule- surface  collisions  are  more  numerous  than  incident 
molecule-incident  molecule  collisions  near  the  body;  and  2.  it  is  assumed 
that  the  incident  and  reflected  molecules  do  not  interact  (collide)  with  each 
other  near  the  body,  i.e.  that  collisions  of  the  incident  molecules  with  the 
surface  are  much  more  numerous  than  incident  molecule-reflected  molecule 
collisions.  Since  the  mean  free  path  by  definition  is  the  average  distance 
travelled  by  a  molecule  between  collisions  with  other  molecules,  a  Knudsen 
number  based  on  Xci  will  indicate  when  the  first  assumption  is  satisfied 
and  a  Knudsen  number  based  on  A^r  or  Aft  will  indicate  when  the  second  as¬ 
sumption  is  satisfied.  To  obtain  a  conservative  estimate  of  when  the  second 
assumption  is  fulfilled,  the  smaller  of  the  two  mean  free  paths,  Af/-  and 
\ri  ,  will  be  used.  It  will  be  shown  that  Ai’r  is  smaller  of  the  two 
mean  free  paths,  and  that,  for  a  near  earth  satellite,  the  Knudsen  number 
based  on  A  if  is  an  order  of  magnitude  smaller  than  the  Knudsen  number 
based  on  A  Li  when  complete  thermal  accommodation  is  assumed.  Thus, 
the  criteria  for  the  existence  of  a  free  molecule  flow  is  that  the  Knudsen  num¬ 
ber  based  on  Afr  rnust  be  greater  than  10. 

The  mean  free  path  is  defined  as  the  average  distance  a  molecule  travels 
in  unit  time  divided  by  the  number  of  collisions  it  undergoes  in  unit  time. 
Relative  to  a  coordinate  system  fixed  in  the  body,  the  preceding  mean  free 
paths  are  given  by  the  following  expressions, 

A  ■  =  %- 

u  &u 

-  1 

Air  “ST 
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°h 

e ri. 

StL. 


where  0j>  is  the  collision  frequency  of  one  type  J  molecule  with  all  the 
type  i 5?  molecules  and  C=  ( J  is  the  average  speed  of  a  molecule. 

At  this  point  it  is  convenient  to  introduce  a  quantity  which  we  will  need  in  the 
future,  X  co  >  t*le  mean  free  path  of  the  incident  molecules  relative  to  them¬ 
selves  referred  to  a  coordinate  system  that  is  fixed  in  the  gas.  This  is  the 
mean  free  path  that  is  tabulated  in  most  atmosphere  tables,  for  example,  the 
1959  ARDC  atmosphere.  From  Eq.  (12A)  Appendix  A,  this  is  given  as 


A  00 


/ 


SPniird' 


CO 


The  collision  frequency  is  independent  of  the  coordinate  system  in  which 
the  motion  of  the  molecules  is  observed.  It  can  be  shown  that  the  collision 
frequency  for  a  gas  with  mass  motion  is  the  same  as  that  for  a  gas  with  no 
mass  motion.  From  eqn.  (10A)  Appendix  A, 

\‘/Z  _ 

(*) 


~  J? (l Jzf  ni  rrd£cc 


a_ 


©rr  —  J~2-  <-7- 


(V 


Unfortunately,  it  is  difficult  to  evaluate  exactly  the  collision  frequency  for  a 
mixture  of  two  gases  which  are  at  different  temperatures  and  one  of  which 
has  a  mass  motion  while  the  other  docs  not.  This  situation  occurs  only  in  the 
vicinity  of  a  body  in  a  free  molecule  flow.  The  collision  frequencies  ©£y 
and  0pc  are  given  approximately  as 


6tr  -  (4) 

Qri  =  TTctrv^  (s) 
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The  tlf  which  appears  in  these  expressions  is  the  true  number  density  of 
reflected  molecules  in  the  vicinity  of  the  body  whereas  the  used  in 

Section  II  B  was  a  fictitious  number  density  which  does  not  exi st  anywhere  in 
the  flow. 

Before  we  can  decide  which  of  the  two  meanfree  paths,  y^ir  or  Art  > 
is  the  smaller,  n^  must  be  computed.  This  will,  in  general,  be  a  function 
of  location  on  the  surface  of  the  body  and  distance  away  from  the  surface  of 
the  body.  Thus,  we  will  compute  an  average  .  From  continuity  considera¬ 
tions,  we  know  that  hT  =  N^.  To  determine  n  ,  we  need  to  know  the  average 
volume  that  the  number  of  molecules  reflected  from  unit  area  per  unit  time 
occupy.  This  is  given  by  an  average  velocity  of  the  reflected  molecules  times 
unit  time,  times  unit  area.  Thus, 

s  _  NrtA.A  _  Nh 

f  VtdA  "  T 

To  see  that  the  total  number  of  molecules  that  wiil  be  contained  in  this  average 
volume  is  N  ,  consider  the  following.  By  the  assumption  of  diffuse  reflection, 
the  same  number  of  molecules  with  a  particular  velocity  will  leave  an  element 
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of  area  in  one  direction  as  leave  it  in  the  opposite  direction.  Thus,  it  would  be 
expected  that  not  all  N?  molecules  would  be  in  the  average  volume  Vt'dA 


surrounded  by  other  elements  of  area  from  which  molecules  are  leaving  in  all 
directions,  it  seems  probable  that  as  many  molecules  will  enter  the  average 
volume  above  dA  in  unit  time  from  the  surrounding  dA's  as  the  number  of 
molecules  which  impinge  on  dA  leave  it  (see  figure  V-Z). 

Since  =  N.,  and  introducing  the  molecular  speed  ratio,  we  obtain 
from  Eq.  (15 B) ,  appendix  B 


There  are  many  average  velocities  that  may  be  defined.  There  are  the 
average  speed  of  a  molecule  which  may  be  moving  in  any  direction  in  space, 
the  average  speed  of  a  molecule  which  can  only  move  in  the  space  above  the 
element  of  area^and  the  average  component  of  velocity  normal  to  the  element 
of  area.  It  turns  out  that  these  velocities  differ  from  each  other  by  less  than 
Z5%.  Since  the  purpose  of  these  considerations  is  to  determine  if  a  free  mole¬ 
cule  flow  exists  and  since  it  is  desirable  to  be  conservative  in  our  estimate, 

we  want  to  choose  the  smallest  reasonable  value  of  which  means  that  n 

At"  r 

should  be  as  large  as  possible,  and  thus  the  average  velocity  we  use  should  be 
the  smallest  one.  The  smallest  average  velocity  is  also  the  one  which  it 
seems  reasonable  to  use,  the  average  velocity  of  the  reflected  molecules 
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normal  to  the  surface.  To  compute  this,  we  proceed  as  follows.  (J~  is  the 
velocity  component  normal  to  the  element  of  area  where  the  positive  normal 
is  directed  inward.  The  number  of  molecules  in  a  given  velocity  range  that 
leave  unit  area  of  the  surface  in  unit  time  is 

(-0-)nr-fa{uct(rJ  Ur- 

where  n^  is  the  number  density  of  the  fictitious  gas  from  which  the  reflected 
molecules  are  assumed  to  issue.  The  velocity  normal  to  the  surface  that  the 
number  of  molecules  in  a  given  velocity  range  that  leave  unit  area  in  unit  time 
carry  with  them  is 


ir(-ir)n^  -fdtu.d.V'd.us' 


The  average  normal  velocity  of  the  reflected  molecules  is  then 


77  _  /  /  / 


V  = 


AJ-  [ZirR.Tr 


Vt 


+  <X>  0  .  ,  7  V  o'! 

r  J—  ( k+Lr-^+Hr7-/ 


dlr-Jur 


- cO-OD-CO 


Performing  the  indicated  integrations  we  find  that 

\7  —  _  ^  2.  RTj~ 

v  ^  ir¬ 


es) 


From  Section  II,  Part  B,  Eq.  (17)  gives  n^  when  the  molecular  speed 
ratio  is  introduced  as 


'V  =  [fiT  )f  s  ( I+Mjis)  +  e  *  s 


(10) 
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Substituting  Eqs.  (8)  and  (10)  into  (9)  gives 

V=  --gfzVRT,.'  (H) 

The  minus  sign  means  that  this  V  is  in  the  negative  y  direction,  i.e.  it  is 
directed  outward,  normal  to  the  surface.  Thus,  we  will  use  only  the  magnitude 
of  V, 

y  =  i  JzVfiTr'  (It) 


Substituting  Eqs.  (8)  and  (12)  into  Eq.  (7)  gives  for  the  average  number 
density  of  reflected  molecules 


nifT 


fjr  rs(if^irs)  + 


Since  we  are  interested  in  the  greatest  density  of  reflected  molecules  which 
can  occur,  we  want  the  maximum  value  of  the  above  expression  for  nf.  This 
occurs  for  large  S  when  — *■  1  which  occurs  in  the  stagnation  region  of  the 

body.  Thus,  in  the  stagnation  region 


ms  nr 

W1  J  Tr 

i  the  expression  fi 

IT 

2  n- 5  Td&0 

'v 

(14) 


oe) 


Forming  the  ratio  and  introducing  the  molecular  speed  ratio  in  A  ri 


gives 


_Ai£  —  If  ^  I 
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Thus,  \ir<\r  .•and  \lr  will  be  used  to  compute  the  Knudsen  number  which 
indicates  when  the  second  assumption  of  free  molecule  flow  is  satisfied.  In¬ 
troduce  A 05  into  Eq.  (i5).  This  gives  for  Av  ^r> 

\  _  -L  fT  i  /  v 

^ it  0  2  sj  -ft  (f*) 

Now  form  the  ratio 

-Air  =  JELflP  in) 

Xu'  sZ  y  77 

From  this  expression,  we  sec  that  for  S  y  1 ,  T  <  T.  (which  is  the  case  for 
a  satellite),  and  nearly  complete  energy  accommodation,  Xi/-*'  At' f  and  the 
critical  Knudsen  number  is  based  on  .  When  this  is  greater  than  ten, 

the  Knudsen  number  based  on  \  ■  £  will  also  be  greater  than  ten.  If  there  is 
an  appreciable  lack  of  energy  accommodation,  tire  critical  Knudsen  number 
may  be  based  on  Xli  •  At  the  time  this  is  written,  it  appears  reasonable 
to  assume  nearly  complete  energy  accommodation  and  to  require  that  the 
Knudsen  number  based  on  Afr  130  Mrttal-cr  than  ten  for  a  free  molecule  flow 
to  exist,  i.  e.  , 

Kn  -  Ah  ?/<? 


For  large  molecular  speed  ratios,  X  it-  becomes  a  severe  criteria  for  de¬ 
termining  whether  or  not  a  free  molecule  flow  exists.  In  marginal  cases,  the 
physical  processes  that  arc  occurring  should  be  given  a  more  detailed  exami¬ 
nation. 

The  preceding  choice  of  \it~  as  the  critical  mean  free  path  differs  from 

Reference  7  in  which  Art  was  chosen.  The  result  of  Reference  7  was  based 

on  the  fact  that  in  this  article  it  was  shown  that  the  actual  number  density  of 

reflected  particles  falls  off  very  rapidly  away  from  the  body,  and  thus,  only 

very  near  the  surface  will  the  number  density  approach  n^  and  \lj~  be  less 

than  Art  •  Since  this  work  was  based  on  the  hyperthermal  assumption  and 

fi  \ «  p  fictitious  density  that  docs  not  exist  anywhere  in  the  flow,  it  wets  felt 
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that  the  use  of  an  average  number  density  as  presented  here  provides  a  more 
appropriate  criteria  for  determining  when  a  free  molecule  flow  exists.  It 
should  be  pointed  out,  however,  that  the  ratio  of  to  \  f-  L  *s  which 


is  almost  one. 


B.  COMPILATION  OF  THE  AVAILABLE  EQUATIONS 


The  equations  for  the  aerodynamic  coefficients  which  are  available  are: 
FLAT  PLATE  WITH  BOTH  SIDES  EXPOSED  TO  THE  FLOW. 


The  orientation  of  this  plate  is  given  in  Figure  III- 1 . 

c<  ■  t, z' 

C„=f  C.  +  fi+tj r„ 


FLAT  PLATE  WITH  ONE  SIDE  EXPOSED  TO  THE  FLOW 
The  orientation  of  this  plate  is  given  in  Figure  III- Z . 

C  -  A  rcir 


c*  =  ^f 


CF  ,\Tr 


X,  + 


r£ 


C3) 


it  snould  be  noted  that  the  center  of  pressure  of  a  fiat  plate  is  always  the 
centroid  of  the  area  of  the  plate.  It  is  usually  more  convenient  to  use  this 
fact  to  compute  the  moment  than  to  use  the  moment  equations  which  are  given 
for  the  flat  plates. 
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RIGHT  CIRCULAR  CYLINDER 

The  orientation  of  the  cylinder  is  given  in  Figure  III- 3. 


FRUSTUM  OF  A  RIGHT  CIRCULAR  CONE 

The  orientation  of  the  cone  is  given  in  Figure  III-4. 
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SPHERICAL  SEGMENT 


The  orientation  of  the  spherical  segment  is  given  in  Figure  III- 5. 

■  2  f  __  5 P  r=r~>  _ sp 


c 


N 


Am/  l 


Ji 

Tt 


The  preceding  equations  are  valid  for  angles  of  attack  which  are  positive 
as  defined  in  the  figures  giving  the  orientation  of  the  different  configurations. 
For  negative  angles  of  attack, 

~CN(oc) 

C/v)(-°0= 

but 


where 


cA(-*)= 


C.  RESOLUTION  OF  A  VEHICLE  INTO  SIMPLE  SHAPES 


The  next  question  which  arises  is  that  of  shielding.  The  equations  which  have 
been  presented  do  not  take  into  account  the  shielding  of  one  part  of  a  vehicle 
by  another  part.  The  only  shielding  they  consider  is  that  of  an  exposed  part 
shielding  a  portion  of  itself  such  as  the  shielding  of  the  back  side  of  a  cylinder 
by  its  front  side.  The  shielding  at  a  given  angle  of  attack  is  best  determined 
graphically  using  the  techniques  of  descriptive  geometry.  For  each  angle  of 
attack,  the  exposed  parts  of  the  vehicle  can  then  be  divided  into  flat  plates, 
cylinders,  cones  and  spherical  segments.  The  aerodynamic  coefficients  are 
computed  at  each  aiigie  of  attack  ior  piece  of  ihe  vehicle  and  are  then 
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added  together  to  give  the  total  coefficients  for  the  vehicle  at  each  angle  of 
attack. 

D.  DETERMINATION  OF  THE  PARAMETERS  S,  T.  AND  T  . 

’  1  r 

The  molecular  speed  ratio  is  given  by  S  -  Vjs*  Ti  where  q  is  the 
vehicle  velocity,  T^  the  temperature  of  the  incident  molecules,  and  R  the  gas 
constant  for  the  particular  gas  of  interest. 

T^  is  the  temperature  of  the  incident  molecules  and  is  obtained  from  an 
atmosphere  table. 

The  value  to  be  used  for  Tf  depends  on  the  degree  of  thermal  accommo¬ 
dation  that  occurs.  The  scarcity  of  experimental  data  for  large  molecular 
speed  ratios  makes  it  difficult  to  specify  the  value  of  the  thermal  accommoda¬ 
tion  coefficient.  In  the  results  presented  in  this  report,  T^  was  assumed  con¬ 
stant  over  the  surface.  Thus,  if  T  actuallv  varies  over  the  surface,  an  aver- 
age  T^  should  be  used.  If  complete  thermal  accommodation  is  assumed, 

Tr  =  Tw  and  the  average  wall  temperature  over  the  surface  should  be  used. 

In  this  case  the  value  of  is  not  very  sensitive  to  T^  because  T^ 

is  of  the  order  of  2200°  R  and  T  is  of  the  order  of  500°R  and  a  ^50°  varia- 

. _ w 

tion  in  T^  does  not  affect  yj  Tr  /Tc  very  much.  Further,  the  value  of  the 

US,  <*  ,  f3  )  due  to  the  reflected  molecules  which  is  multiplied  byj7 ?/t/ 

is  usually  an  order  of  magnitude  smaller  than  thel(  S  ,  0(  ,  f&  )  due  to  the 
incident  molecules.  Thus,  for  complete  thermal  accommodation,  the  error 
introduced  into  the  aerodynamic  coefficients  by  using  an  average  Tw  is  not 
very  large.  At  the  present  time,  it  appears  reasonable  to  assume  nearly  com¬ 
plete  energy  accommodation. 

If,  however,  very  little  thermal  accommodation  occurs,  the  value  of  T^ 
may  be  very  large  due  to  the  large  mass  velocity  of  the  incident  molecules 
(of  the  order  of  25,  000  ft/ sec)  which  gives  them  a  large  kinetic,  energy.  Then, 
if  there  is  little  thermal  accommodation,  the  energy  of  the  reflected  mole¬ 
cules,  which  for  diffuse  reflection  is  measured  by  T^,  will  be  very  nearly 
equal  to  the  energy  of  the  incident  molecules  and  thus  T  will  be  large.  In 
fact,  under  these  conditions  T^  T  and  the  value  of  T^  will  significantly 
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affect  the  value  of  the  aerodynamic  coefficients.  In  this  case,  the  value  of  T 
must  be  obtained  from  an  energy  balance  for  the  surface.  Thus,  it  is  impor¬ 
tant  that  the  latest  results  of  work  on  the  problem  of  thermal  accommodation 
be  incorporated  as  soon  as  possible  into  the  free  molecule  aerodynamic  analy¬ 
sis  . 


E.  PITCH  AND  YAW  COMPONENTS  OF  THE  AERODYNAMIC 
COEFFICIENTS 

The  equations  presented  here  give  the  aerodynamic  coefficients  in  the 
plane  containing  the  velocity  vector  and  the  longitudinal  axis  of  the  vehicle, 
the  moment  being  about  the  axis  normal  to  this  plane.  Thus,  the  forces  in  and 
the  moment  normal  to  the  plane  containing  the  velocity  vector  and  the  longitu¬ 
dinal  axis  can  be  determined  for  an  arbitrarily  oriented  vehicle  in  a  uniform 
free  molecule  stream.  It  is  sometimes  convenient  to  consider  the  pitch  and 
yaw  components  of  the  forces  and  moments  on  a  vehicle.  Since  the  equation 
for  the  force  on  an  element  of  area  is  non  linear,  this  situation  must  be  handled 
by  computing  the  forces  in  and  moment  about  an  axis  normal  to  the  plane 
formed  by  the  mass  velocity  vector  and  the  longitudinal  axis  of  the  vehicle. 

The  and  are  then  resolved  into  a  pitch  component  and  a  yaw  component, 
is  the  total  and  remains  unchanged. 

The  values  of  the  free  molecule  flow  functions,  I  (  S  ,  Of  ,  {3  )  are  ob¬ 
tained  from  the  tables  once  the  molecular  speed  ratio  S,  angle  of  attack  CX  , 
and,  in  the  case  of  the  cone  and  spherical  segment,  the  semi  vertex  angle  6 
or  the  sphere  half  angle  ©  ^  are  known. 


A 

IS 
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SECTION  VI 

CONCLUDING  REMARKS 

At  the  present  time  there  has  not  been  an  experimental  determination  of 
the  free  molecule  aerodynamic  coefficients  of  any  bodies  for  Knudsen 
numbers  greater  than  10  and  molecular  speed  ratios  between  5  and  10.  The 
only  experimental  work  which  has  been  done  in  the  free  molecule  regime  is 
the  heat  transfer  and  drag  measurements  of  Reference  (8)  for  a  cylinder 
transverse  to  the  flow.  In  this  experiment  the  Knudsen  number  varied  from 
4  to  1 8 F  and  the  molecular  speed  ratio  varied  from  0.5  to  2,3.  The  high 
Knudsen  numbers  were  obtained  by  placing  a  0,0031  inch  diameter  wire  in  a 
low  density  wind  tunnel. 

The  lack  of  experimental  data  is  due  primarily  to  two  things.  First, 
it  is  very  difficult  to  obtain  a  free  molecule  flow  for  which  the  Knudsen 
number  is  greater  than  10  and  still  use  a  model  which  is  large  enough  to  be 
used  for  force  and  moment  measurements,  and  second,  even  if  a  large 
model  can  be  used,  the  magnitude  of  the  forces  to  be  measured  is  within  the 
limits  of  accuracy  of  available  equipment.  Thus,  there  is  essentially  no 
experimental  data  available  that  can  be  used  to  check  the  theoretical  pre¬ 
dictions  in  the  Knudsen  number  and  molecular  speed  ratio  range  of  interest. 

To  further  complicate  the  picture,  tne  molecule -surface  interaction 
phenomena  and  its  variation  with  angle  of  incidence,  molecular  speed  ratio, 
temperature,  etc.  is  not  very  well  known.  There  is  also  the  problem  of  the 
shielding  of  one  surface  by  another,  and  the  reflection  of  molecules  from  one 
surface  onto  another.  This  latter  problem  has  been  discussed  and  results 
given  for  the  sphere  (references  9  &  10)  and  infinite  cylinder,  reference  9. 

It  is  evident  that  there  is  a  large  amount  of  work  that  must  be  done  be¬ 
fore  we  will  have  a  complete  understanding  of  the  processes  occurring  in  a 
free  molecule  flow.  For  these  reasons,  it  is  not  possible  at  this  time  to 
estimate  the  accuracy  of  the  theoretical  free  molecule  aerodynamic  coef¬ 
ficient  predictions.  The  only  thing  that  can  be  said  is  that  they  a1*.-  accu¬ 
rate  as  we  can  predict  theoretically . 
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APPENDIX  A 

FUNDAMENTAL  ASSUMPTIONS  AND  IDEAS 
OF  THE  KINETIC  THEORY  OF  GASES 

The  kinetic  theory  of  gases  rests  upon  two  closely  related  hypotheses: 
The  first  is  that  of  the  molecular  structure  of  matter  and  the  second  is  that 
the  degree  of  molecular  motion  is  determined  by  the  heat  content  of  the 
body,  that  is,  the  internal  energy  of  the  gas  is  contained  in  the  random 
translational  motion  of  the  molecules.  The  physical  model  upon  which 
kinetic  theory  is  based  is  that  of  an  assembly  of  gas  molecules  in  a  con¬ 
tainer  in  random  motion,  bumping  into  themselves  and  the  walls  of  the  con¬ 
tainer.  Kinetic  theory  is  a  statistical  theory,  that  is,  it  is  based  on  the 
laws  of  chance . 

Theoretically  it  is  possible  from  the  laws  of  dynamics  to  completely 
describe  the  motion  of  all  the  molecules  of  a  gas  in  a  container.  However, 
this  is  of  little  interest.  Complete  knowledge  of  the  history  of  every 
molecule  in  the  gas  is  not  what  we  are  seeking;  we  need  to  know  only  so 
much  about  the  molecular  motions  as  will  suffice  to  explain  and  predict 
quantitatively  the  observable  macroscopic  properties  of  the  gas  such  as 
pressure  or  viscosity.  Kinetic  theory  attempts  to  explain  these  properties 
in  terms  of  the  molecular  structure  and  the  forces  between  the  molecules. 

To  say  anything  quantitative,  we  must  say  something  about  the  geo¬ 
metry  of  the  molecules.  We  will  picture  a  molecule  as  a  rigid,  perfectly 
elastic,  spherical  body,  that  is,  a  billiard  ball.  This  assumption  means 
that  no  momentum  or  energy  will  be  transferred  to  the  internal  degrees  of 
freedom  of  the  molecule  and  thus  momentum  and  energy  will  be  conserved 
during  collision.  The  justification  for  this  assumption  of  a  billiard  ball 
molecule  is  in  its  success  in  predicting  the  properties  of  a  teal  gas.  One 
of  the  most  striking  features  of  the  kinetic  theory  is  the  extent  to  which  it  is 
possible  to  predict  the  behavior  of  a  gas  as  a  whole,  while  remaining  in 
almost  complete  ign-  ranee  of  the  behavior  and  properties  of  the  molecules 


A-l 


LOCKHEED  AIRCRAFT  CORPORATION 


MISSILES  and  SPACE  DIVISION 


LMSC-448514 


of  which  it  is  composed.  The  reason  for  this  is  that  many  properties  of 
gases  depend  only  on  general  dynamical  principles  such  as  the  conservation 
of  energy  and  momentum. 

Four  things  are  required  to  specify  the  properties  of  a  gas: 

1.  The  diameter  of  each  molecule,  2.  the  mass  of  each  molecule, 

3.  the  number  of  molecules  per  unit  volume  and  4.  the  distribution  of  the 
random  motion  of  the  molecules.  Once  the  particular  gas  of  interest  is 
specified,  the  molecular  diameter  and  mass  are  known.  The  number  of 
molecules  per  unit  volume  is  determined  once  the  density  is  specified. 

That  leaves  only  the  distribution  of  the  random  motion  of  the  molecules  to 
be  determined. 

A  distribution  function  may  be  defined  as  follows  .  A  distribution 
function  gives  some  sort  of  value  (i.e.  density,  velocity,  etc.)  associated 
with  a  "volume"  in  some  "space".  We  want  to  find  the  velocity  distribution 
function  for  a  steady  state  distribution  of  velocities.  In  the  steady  state, 
the  distribution  function  will  be  independent  of  time. 

A  molecule  may  be  said  to  have  six  coordinates,  three  spacial  coor¬ 
dinates,  and  three  velocily  coordinates.  Let  the  spacial  coordinates  be 
x,  y,  z  and  the  velocity  coordinates  be  u,  v,  w  corresponding  to  velocities 
in  the  x,  y,  and  z  directions  respectively.  Let  n  be  the  number  of  molecules 
per  unit  volume.  We  want  to  know  how  many  molecules  have  velocities 
which  lie  between  u  and  u  -f  du,  v  and  v  -f  dv,  w  and  w  +  dw.  If  we  define 
the  density  of  velocities  with  coordinates  u,  v,  w  as  'C'  ,  then  /T*  (u,  v,  w) 
is  the  distribution  function  for  velocity  in  the  velocity  space.  It  gives  the 
number  of  points  in  an  element  of  volume,  du  dv  dw,  in  velocity  space. 

Since  each  point  corresponds  to  one  molecule  in  x,  y,  z  space,  gives 
the  number  of  molecules  with  velocities  in  du  dv  dw.  Then  the  number  of 
molecules  per  unit  volume  with  velocities  between  u  and  u  f-  du,  v  and 
v  +  dv,  w  and  w  +  dw  is 

dn  ~  T' cLu.  dtr dor 

We  will  define  the  normalized  distribution  function  as 
a . \  _  'C  (u>  v,  w) 

4u,  V  ,  v.  J  -  - — - 
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Using  this  notation,  dn  becomes 

cin  =  nf(u-jv-jus-)  ditdircLur 


Since  n  is  the  number  of  molecules  per  unit  volume,  the  probability  that  the 
velocity  of  a  molecule  selected  at  random  from  the  element  of  volume  dx  dy 
dz  will  be  in  the  range  u  to  u-f-du,  v  to  v-#-dv  and  w  to  w*f  dw  is 


dn 

n 


-  f  {u,  v,  v/}  du  dv  dw. 


We  now  search  for  a  steady  state  law  for  the  distribution  of  velocity  in  a 
gas.  A  rigorous  derivation  of  the  Maxwell  velocity  distribution  function  is  a 
tedious  procedure.  Since  its  existence  has  been  verified  experimentally,  it 
may  be  considered  a  known  physical  law  and  the  following  discussion  may  be 
considered  as  a  plausibility  argument  for  its  existence.  Since  it  is  assumed 
that  the  molecules  in  a  gas  do  not  effect  each  other  except  during  collision, 
the  only  way  the  velocity  of  a  molecule  can  be  changed  is  by  a  collision  with 
another  molecule.  Consider  two  molecules  in  collision,  molecules  A  and  B 
with  velocities  before  and  after  collision  as  indicated  below. 


Molecule  A 
Molecule  B 


Before  Collision 


After  Collision 


u,  v,  w 

t  f  r 
U  ,  V  ,  w 


U,  V,  w 

_ f  _ / 
U  ,  V  ,  w 


In  order  to  determine  the  velocities  after  collision,  the  direction  of  the  line 
of  centers  as  well  as  the  initial  velocities  must  be  known.  Let  type  (X 
collisions  be  collisions  of  molecules  with  given  velocities  (u,  v,  w)  and  (u  , 
v/,  w f )  and  given  direction  of  line  of  centers.  The  number  of  collisions  of 
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this  type  will  be  proportional  to  the  number  of  molecules  of  class  (u,  v,  w) 
and  of  class  (u^,  v  /t  w^)  which  are  present. 

The  number  of  molecules  of  class  (u,  v,  w)7  -  nf  (u,  v,  w)dudvdw 

per  unit  volume.  J 

The  number  of  molecules  of  class  (u^,  v/,  w/)L  =  nf  ( u  ^  v/>  w;  ')du'<l<rkur 
per  unit  volume.  J 

The  number  of  collisions  of  type  Q(.  will  be  proportional  to  the  product  of  the 
number  of  molecules  of  classes  A  and  B  that  are  present.  Thus  the  number 
of  collisions  of  type  Q(  per  unit  volume  per  unit  time  is 

where  A^  is  the  function  of  proportionality.  A^  is  a  function  of: 

(1)  the  diameter  of  the  molecules  -  the  larger  the  molecules, 
the  more  likely  they  are  to  have  a  collision. 

(2)  the  relative  velocity  of  the  molecules  -  molecules  moving 
parallel  to  each  other  are  not  likely  to  run  into  each  other 
while  if  they  are  moving  toward  each  other,  they  are  very 
likely  to  have  a  collision. 

(3)  the  direction  of  the  line  of  centers. 

To  be  completely  rigorous,  we  should  determine  the  functional  form  of  A^  . 
However,  this  is  a  detailed  question  and  we  will  not  go  into  it  since  it  is  not 
necessary  to  know  A  ^  in  order  to  determine  the  distribution  function. 

During  each  collision  of  type  CX  ,  one  molecule  of  class  (u,  v,  w)  and  (u/, 
v w/)  is  lost  and  one  molecule  of  class  (u,  v,  w)  and  (u/,  v  ,  v//)  is  gained. 

Now  let  a  type  {3  collision  be  the  inverse  of  a  type  Oi,  collision,  that 
is,  it  is  a  collision  with  velocity  (u,  v,  w)  and  (u  ,  v  ,  w  )  before  collision 
and  velocity  (u,  v,  w)  and  (u/,  v  wO  after  collision.  Using  reasoning 
similar  to  that  for  a  type  (X.  collision,  we  find  that  the  number  of  collisions 
of  type  yS  per  unit  volume  per  unit  time  =  n  f  (u,  v,  w)f(u  ,  v  ,  w^) 

dudvdwdu  ,  dv  '  dw  /  . 

The  relative  velocities  of  the  molecules  in.  type  (X  and  (3  collisions 
are  the  same.  The  diameters  of  the  molecules  are  the  same  because  all  the 
*  molecules  are  physically  the  same.  From  the  dynamics  of  the  collisions  and 
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the  definition  of  a  type  y?  collision,  the  direction  of  the  line  of  centers  is  the 
same  for  type  p(  and  [3  collisions.  Therefore, 


A  c<  ~  A  (i  , 

It  can  be  shown  that  dudvdw  du',dv/dw'’  =  dudvdw  du  dv  dw  .  Using 
these  facts, 


the  number  of  coll 
per  unit  volume  pe 


isions  of  type/3]  _  A  r\Z-f(Uj  ir)  Ur)f(£j  ir} w~' )dudtr±wju  'Jtr'eJur 
:r  unit  time  J  ^ 


The  net  gain  of  molecules  of  cla  ss  (u,  v,  w)  per  unit  volume  per  unit 
time  as  a  result  of  collisions  of  types  pi  and  /3  is  equal  to  the  number  of 
molecules  gained  by  type  [3  collisions  minus  the  number  of  molecules  lost 
by  type  collisions 

=  A^  n^£ffu,  v,  w)f(u  ,  v  ,  w  )-f(u,  v,  w)f(u/,  v'f  w/)jdudvdw  du**  dv/  dwy 

If  we  set  the  term  in  brackets  equal  to  0,  we  have  that 

f(u,  v,  %v)f(u  ,  v  ,  w  )-f(u,  v,  wjfti/,  v',  w' )  -  0.  (2-A) 


This  means  that  there  will  be  no  net  gain  or  loss  of  molecules  of  class  (u, 
v,  w).  If  this  term  is  zero,  then  since  A  ^  which  depends  on  the  direction 
of  the  line  of  centers  does  not  appear  in  the  brackets,  it  is  zero  regardless 
of  the  direction  cf  the  line  of  centers.  Further,  if  we  require  that  the  term 
in  brackets  be  zero  for  all  velocities,  then  the  gas  is  in  equilibrium. 
Equation  (2 -A)  is  a  sufficient  condition  for  equilibrium.  Boltzmann  showed 
that  it  is  also  a  necessary  condition  for  equilibrium.  Thus,  for  a  steady 
state  to  exist,  every  collision  and  its  inverse  must  be  in  equilibrium.  This 
is  a  special  case  of  the  principle  of  detailed  balance  which  states  that  if  a 
substance  is  in  macroscopic  equilibrium,  every  microscopic  process  occur- 
ing  inside  the  medium  must  be  in  equilibrium. 

In  order  to  determine  the  equilibrium  velocity  distribution  function,  we 
must  solve  Equation  (2-A)  subject  to  the  restraints  placed  on  u,  v,  w,  u' , 
vJ  ,  w/,  u,  v,  w,  u^,  v/,  w  by  the  energy  and  momentum  conservation 
equations  for  a  molecular  collision.  Take  logarithms  of  Equation  (2A)and 
rearrange.  This  gives 

1  n  f(  u ,  v ,  v.' )  In  f|u^,  v  ‘ ,  w '  )  =  1  n  f(  u ,  v ,  w)  +  In  f(u  ,  v  ,  w  ) 


A-  5 


LOCKHEED  AIRCRAFT  CORPORATION 


MISSILES  and  SPACE  DIVISION 


LMSC-448514 


Let  F  =  Inf.  Then,  F(u,  v,  w)«fr,F(u/,  v;  ,  w/)=F(u,  v,  w)  '4"F(u  ,  v  ,  w  ) 

(3-A) 

To  solve  equation  (3A)  we  need  some  functions  of  velocity  which  have  the 

same  form  before  and  after  collision.  We  have  four  such  functions,  three 

from  the  conservation  of  momentum  and  one  from  the  conservation  of  energy 

2  2  2 

(mu,  mv,  mw  and  m(u“-j-  v  ~/-w  ).  Since  each  of  those  is  a  solution  of  the 
equation,  a  linear  combination  of  them  will  also  be  a  solution.  Thus 

Z  2  2 

F(u,  v,  w)  =  <y,  m(  u  -f  v  -f  w  )  +  CX  ^  mu  mv  -f*  (X  ^  mw  ~h  . 

It  can  be  shown  that  this  is  the  only  solution  of  this  equation.  Rearranging 
the  constants  we  obtain 

In  f(u,  v  w)  -  CX  ,m  f  (u-u  )2-f-  (v-v  )2-l-  (w-w  1  +  6 

’  t  [_  o  o  o  J 

or 

f  (u,  v,  w)  '=  Ae“i®  ^  [  (u-uq)2  4*  (v -vq)2  (w-wq)2  J  (4 -A) 

We  now  have  to  determine  the  constants  which  appear  in  equation  (4-A), 
The  probability  of  finding  a  molecule  for  which  (u-uq)  is  negative  is  the  same 
as  the  probability  of  finding  a  molecule  for  which  (u-uq)  is  positive,  Thus, 
the  average  value  of  (u-uq)  is  zero.  Since 

u-u  =  0,  u-u  =  0 

o  ’  o 

and  u  =  u  =  u  because  u  is  a  constant.  This  indicates  that  u  is  the  mass 
o  o  o  o 

velocity  of  the  gas.  Thus,  u  ,  v  and  w  are  the  components  of  the  mass 

’  o’  o  o 

motion  of  the  gas.  For  no  mass  motion, 

-  /3m(u>f  v^fw2) 

f(u,  v,  w)  =  Ae 

where  in  this  case,  u,  v  and  w  are  the  components  of  the  thermal  motion  of 
the  molecule  s  , 

Now,  the  number  of  molecules  with  a  particular  velocity  is  d.V~dur. 

The  integral  of  this  over  all  possible  velocities  will  give  the  total  number  of 
molecules  in  a  unit  volume.  Thus 
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+  00 


jj' f  dudvdw  =  n 


—bO 

Since  n  is  a  constant,  it  can  be  factored  out  of  the  integral.  Then 

+  00 


F  dudvdw  =  1 


-  OO 


Substituting  in  for  f  in  the  case  of  no  mass  motion  gives 


+00  +  co  +00 

Afe-^Je^jJe-^dvr  =, 

-00  -  cO 


Performing  the  integrations  we  obtain 


It  can  be  shown  that  the  average  value  of  u“  is 

li  -  = 


where  k  is  the  universal  gas  constanL  per  molecule. 
The  average  value  of  u^  can  be  determined  by 

+  0 0 


ur)du.d(+dur 


~  00 


(S-A) 
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Equating  these  two  expressions  for  u  gives 


+O0 


+9* 


tf  =  Ar  =  /x 


e.  dir 


-mfiitr 

e  ctw~ 


-  00 


_l _ 


/ 

•  r  '  ^KT 


(&-a) 


Substituting  this  value  for  into  the  expression  for  A  gives 

-3/Z 


[i-fC] 


Using  these  values  for 


P 


-f-iUjlTjUr)- 


(JH - 

[t'rrKT 


and  A,  f  becomes 

3/2  _  an —  /  £4e  •+  ir  ^~h  u/~^ 

I  e  SKT 


(*-A) 


l 

i 


This  is  the  Maxwell  velocity  distribution  function  for  a  gas  in  equilibrium 
■with  itself  where  u,  v,  and  w  are  the  components  of  the  random  thermal 
motion  of  the  molecules  , 

Since  we  must  know  the  mean  free  path  of  the  molecules  in  order  to 
determine  the  flow  regime,  we  will  now  derive  an  expression  for  this  quan¬ 
tity,  Consid  r  the  frequency  of  collisions  be'-vv.e"  melrculee  of  type  1  and  2. 
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We  have  a  volume  V  with  molecules  of  type  1  with  velocity  denoted  by 
(uj,  Vj,  w^)  and  molecules  of  type  2  with  velocity  denoted  by 
C2  ~  C2  ^U2’  V2’  w2^  £recluency  o£  collisions  will  depend  on  the  relative 

velocity  between  molecules  1  and  2,  If  the  relative  velocity  is  zero,  the 
molecules  will  not  collide.  The  relative  velocity  is  defined  as 


c 

r 


u2)2  -f-  (vj 


v?  )2  "f* 


1 


We  want  the  average  c^  over  all  molecules  of  typejl  and  2.  First  we  will 
average  c^  for  a  given  type  1  molecule  over  all  type  2  molecules.  Then  we 
will  average  this  over  all  type  1  molecules  to  obtain  the  average  over  all 
molecules  of  types  1  and  2, 

The  number  of  molecules  of  type  2  with  a  given  velocity  is 


N2  *C2^  fiu2f)v2flw2* 


The  average  speed  of  one  molecule  of  type  1  with  velocity  c^ 

molecules  of  type  2  is 

+  00 


-■AS- - 


N2  f2(c2)  du2dv2dw2 


N, 


relative  to  all 


The  number  of  molecules  of  type  1  !■  Z  Nj  f^  (u^,  v^,  w^)  dujdv^dWj, 
with  given  velocity  Cj  J 

}* 


The  average  velocity  of  all  molecules^  _  — 
of  type  1  relative  to  all  molecules 
of  type  2 
-f-cO 
'  '  f r 


JJ 


c  fn  (u,v,wJ  du,dv.,dw_. 
r  Lull  L  L  L 


N  j  fj  (uj  Vj  Wj)  du^Vjdwj 


c  —  / 

r  J 


•+■  00 


N, 


c  f.  (u.v  w  )  f  ju.v  w  )du,dv.dw,du,dv,dw- 
r  1  1  1  -I  1 j  n  C  Z  1  1  ILL  U 


—  CO 
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It  can  be  shown  that  in  a  mixture  of  gases  in  equilibrium,  each  gas  has  a 
maxwellian  velocity  distribution.  Thus 


c,  = 


3/£ 

(hr),  mz) 

(Z  Trkrf 


-too 


~  ZfCT 


(m,  cf-tr^z  cz ) 

d.U,d<rld^dU-zdtridlc(rt 


—  oo 


Performing  the  indicated  integrations  we  obtain, 


8  KT  \ 
irmZj 


(3-A) 


x 

where 


m,  mz 
m,  +  mz 


is  the  reduced  mass. 


We  now  wish  to  find  the  frequency  of  collision  between  type  1  and  type  2 
molecules.  Consider  a  type  1  molecule  moving  through  type  2  molecules. 


Figure  A-2. 


The  radius  of  the  sphere  of  influence 

__  <L ,  +dz  , 


is  d 


12 


The  type  1  molecule  is  moving 
relative  to  the  type  2  molecules  with 
ihe  average  relative  velocity  c  . 


Crt 


/ 

/ 

T~ 

\ 


\ 


1 

/ 


Figure  A-3 

The  distance  traveled  by  a  type  1  molecule  in  time  t  is  c^t. 


The  volume  swept  out  by  the  sphere  of  influence 
per  unit  time 


]•  =  Ch7Td 


i 

/2 
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The  average  rate  of  collision  of  one  type  1  molecule  with  type  2  molecules 
is  Then 


©,*  ®  *rd\  Cr  «  rn&  d; 


'%KT  \^ 

rn*) 


This  holds  for  any  number  of  molecules  1,  2,  3, . S  , 

Thus,  the  average  rate  of  collision  of  one  type  1  molecule  with  molecules 
of  types  1,  2,  3 . Sis 


©<  =  f©,s  =(BVHT)^± 

^  s 


The  mean  free  path  of  type  1  molecules  is  the  average  distance  traveled 
between  collisions.  It  is  given  by  the  ratio  of  the  distance  traveled  per  unit 
time  to  the  average  number  of  collisions  per  unit  time.  Thus^ 


where  c,  is  the  average  velocity  of  the  type  1  molecules. 


'*•  /—  -~r--i  ( H- A) 

*  £  n*  d‘*  J1*  %  ’  ' 

This  expression  is  rigorously  accurate  for  a  mixture  uf  molecules  in  thermo¬ 
dynamic  equilibrium.  If  all  the  molecules  are  the  same,  \  becomes 


X  = 


jtrhd' 
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APPENDIX  B 

integration  of  the  equation  for  the  forcf 

ON  AN  ELEMENT  OF  AREA 


We  want  to  evaluate  the  integrals  in  the 

+.,50  tea  +et>  B  me 

—  =r  m 

dA 


equation 


'ni  f  (ktL+Jtr  -hturj  Lri-dcLcLirdur 

J  J  J 


-CO  0  -CO 


+  mnr 


+•00  0  -too 

r 


j  f  (ka.  +  4(r+t<*r)(r’fdtl<l(rcL 


Ur- 


-oo  -oo 


(/-  S ) 


where 


-  ^r[(^~2j2+((r-LrJ  -{-(ur-arf'J 


f  ~  \SVrt)  e 


(JL)*/2  +(ur~^)Z  ] 


l  W 


c*-®; 


and 


f 

2/?T 


*  *»VV 


, ,  6  - 1  UUL  Wlin  aid  of  appendix  one  i 

ence  (  4).  Consider  the  first  term  of  Eq.  ( I  -  B) .  It  can  be  written  as 
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We  will  evaluate  the  integrals  of  Eq.  (3-B)  separately 


■pac  fpt  +oe 


+•(*  +r>  3/2.  r  r  r 

J"  J UtrfdudtrJuj-  =  (J?pj  J  J  JLUr 

-A  o  -o o  - 


-p[(lt-tL)'L-+(<r-ir)  +(ur-ur)  J 
‘e  dudirdur  (4~B) 


a  -so 


Let  the  components  of  the  thermal  motions  of  the  molecules  be  repre- 


sented  by 

Then 

in  x 

direction, 

H2 

in  the  y  direction  and  in  the  z  direction. 

u 

*  •«  a 

+ 

H1  > 

u  - 

u  =  Hl 

V 

-  V 

+ 

h2  , 

V  - 

v  =  H2 

w 

=  w 

+ 

H3  , 

w  - 

w  = 

du  =  d.H  ^ 
dv  =  dH2 
dw  -  dH^ 
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because  u,  v,  and  w,  the  components  of  the  mass  velocity,  are  constants. 
Making  this  substitution  eq.  (4-B)  becomes 

+*-/-:*> -r*>  3.+*>  V°  t  z  uV 

J  J  !  utr'fdud‘r,iur  =  (y)  j  J  "'IKrt’t///, 


-JO  -~r 


Since  the  random  motions  are  independent  of  each  other,  the  above  equa- 


A  •>  rt* 


tion  can  be  written  as  .  r-+c9  ->  -t-cQ  . 

J  f  =  fifty 


f.  T  -i-aO  _ 

il  ju(r^aciir<^ur  ~  O  +  JfiT 

-cO  O  ~*> 


-(SC, 

~£ 


The  second  integral  is  evaluated  as  follows 


1 


h  .  .  v2  >*>  t  +r°°  C-/3HT 

f  jj  tr 2  fueler  d«s=(£)  (r'+Zirtk+MtJe  ^zje  <*H 3 


-ot>  0  -jo 


=  ($  'ZJf  (f  f ,+ 'V1#,)  f ' 


_  -ys^2 


p»6  +c«  4ot> 


j  J  J+£('W50V  «-•> 


-CO  o  —jo 
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J  p^en'd 

vo  a  ~cA  -**  _p.  —  <*> 

fA  +vO  +oC  _  __  2{jL.Z 

J  ^‘Cr-fdud'-reJcjj-  =  +  jj§  ^  J  (t~A) 


The  third  integral  is  evaluated  as  follows 


0°  _  ,,2. 


-c*  O  ~?Q 


ha  H *  -foO 


-M ] 

2  dH3 


-fO  o  — oo 


Substituting  into  Eq.  (3-B),  equations  S-B),  (6-B)  and  (7-B)  we  get 
■f-A  +tX>  +<*> 


rw^j  j  J (ku  -j-Jtr  +t ur) trfdudirdur 


-0®  O  .00 


=  mnp2jr 


-/3Lre- 

(ku+Jf  +  t&j[lFjT}( !+*^irj3')HjL;e  J 


+  ~  (/+zywj?) 
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Now  consider  the  second  term  of  Eq.  (1-B). 
+  O  ■1L'» 


rY}r>f'  I  (ku-+£u-  -i-tttr)  Lr--fcJuchrc/ur 


It  can  be  written  as 


This  term  gives  the  force  due  to  the  re- emitted  molecules  which  we 
assumed  reflected  diffusely,  i.  e.  they  are  assumed  to  issue  randomly  from  a 
fictitious  gas  which  has  a  maxweilian  velocity  distribution  corresponding  to 
and  no  mass  motion  relative  to  the  element  of  area.  Thus,  the  only  ve¬ 
locity  the  reflected  molecules  possess  is  their  random  thermal  motion  cor¬ 
responding  to  T^.  Thus, 

u  =  Hj,  v  =  w  =  and 

du  =  d.Hp  (iv  =  dH7,  dw  =  dHy 


+00  0  +P° 

111  uu~^ac^^ur 
-<*)  -O0  -00 


-0$ 


-fir"? 


<5 -  X  +*> 


We  will  evaluate  these  integrals  one  at  a  time 

-  ^  f*  r  -Wi, 
dH,\Hze 


=0  (to- 


~oC 


-0 C 


3  j  - 

j  j  j  irZ  fdudir  ctur-  = 


~ab  46  -cb 


O  7  _ 


-od 
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Let  =  -  x  Hz  -  -rfj  X  =  <*> 

dHz  ~  -d)(  tf2=Oj  X  =  O 

^  7  0  p 

,,L  ~  dr  J  f  2_  .  f 

HZC  riH2  ~  ~~J  X  t  cix  =  J 


0°  z 

-Pr* 


p  -1-'  -' ,  /  /  7r 

X  e  ^x=  4/W^r 


-*0 


p6 


+  \>  0 

ir^fjcidirjur  -  r^- 

4-Pr 


~cxS  -<&  -*> 
+ot>  o  +•<> 


~H>C  ■j  ®  . ,  <-  /• 

7  r°  M 


(ii- a) 


-h*>  ? 


=  O 


0*-B) 


Substituting  eqs.  (10-B),  ( 1  1  -  R)  and  (12-B)  into  eq.  (9-B)  gives 

f.<®  Cl  tco 

/77/vJ  j  j  (ktl+Jtr+tur)irfJudirdur  =  Qi-g) 

-io  ~c0  ^00 


where 


Z3^ 


/ 


B-6 


LOCKHEED  AIRCRAFT  CORPORATION 


MISSILES  and  SPACE  DIVISION 


JL.MSC-4485 1 4 


At  first  it  might  seem  surprising  lhat  the  reflected  molecules  contribute 
only  a  force  normal  to  che  surface.  However,  if  it  is  remembered  that  when 
we  assumed  that  the  molecules  where  diffusely  reflected,  it  was  assumed  that 
they  were  reemitted  in  a  random  manner,  that  is,  as  many  molecules  with  a 
particular  velocity  leave  the  surface  in  one  direction  as  in  the  opposite  direc¬ 
tion.  Thus,  the  tangential  momentum  components  cancel  each  other  out. 
Therefore,  this  result  is  to  be  expected. 

Substituting  equations  (8-B)  and  (13-B)  into  equation  ( 1  - B )  gives  for  the 
component  of  the  total  force  on  a  element  of  area  in  a  particular  direction 


-  mnL  p—  (kZ+Jir  +tay-)^fir' (i+W^j+jr  ] 
L  2fTr  u 


dA 


0  4--B) 


Determination  of  n  : 

r 


The  number  of  molecules  in  the  velocity  range  u  to  u  +  du,  v  to  v  +  dv,  w  to 
w  +  dw  that  impinge  on  the  element  of  area  in  unit  time  is 


til  irfciudcrduf 


The  total  number  of  mc.ecules  in  all  velocity  ranges  that  impinge  on  the  ele¬ 
ment  of  area  in  unit  time  is 

j-CO  -hOO  +■*> 

1  /  /" 


A /<;  =  /?’J  J  j  Lrf-dudv-dw- 


-ot  £  -c6 
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This  integral  is  evaluated  in  a  manner  similar  to  the  previous  integrals. 

-f-o o  z  +2*  z  wt 

Nc  -  nL^ j**J £‘3H‘AH,j(o:  + Hz) e  ^ 


/  a  \  VZi — i  r  I — I 

=  rjlf£)  /Ffe  IF 
W  J  P  I'zjp 


(tr-e) 


We  proceed  in  a  similar  manner  for  the  reflected  molecules. 


The  number  of  molecules  in  the  velocity  range  u  to  u  +  du,  v  to  v  +  dv,  w  to 
w  +  dw  that  leave  the  element  of  area  per  unit  time  is 


nh  C-  U-)  fdu  Aird  ur- 


The  total  number  of  molecules  in  all  velocity  ranges  that  leave  the  ele¬ 
ment  of  area  in  unit  time  is 


Nr  =  nr 


f  I  f  (-ir)'f  AilAtrAur 


-To  -oO  -00 


Since  the  fictitious  gas  from  which  the  reflected  molecules  are  assumed 
to  issue  has  no  mass  motion  relative  to  the  surface,  the  integral  becomes 


V*  *v//3 
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L£  r 

r° 

-hf,e 


~cfX 


I  Ht,el’rHljHt  =  -/; 

~C0  00 


dnd. 


Hz  =  -c£> )  X- ob  ■  Hz~  °j  X  =  0 


-fir*, 

X  £  dX  — 


xe  dx  - 


^^3r 


ZjWr 


0<.-e) 


From  continuity,  since  molecules  cannot  be  created  or  destroyed, 


Ni 


(11-0) 


Substituting  in  eq.  (17-B)  for  N.  and  and  solving  for  n^  gives 

nh  =  nL  ^[Trir(i+Mt?J?)+'j=re/3  ] 


(n-B) 


The  function,  erf  x,  which  appears  in  the  preceding  equations  is  defined 


as 


?  f  *  -y2 
e  d 1 


(13-  B) 
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APPENDIX  C 

TABLES  OF  THE  FREE  MOLECULE  FLOW  FUNCTIONS 

Due  to  their  bulk,  the  tables  of  the  free  molecule  flow  functions  have 
been  published  under  separate  cover  and  will  be.  supplied  to  interested  persons 
upon  request.  Requests  should  be  sent  to: 

Project  Publications  and  Services 
Department  61-62 

Lockheed  Missiles  and  Space  Company 
Sunnyvale,  California 

The  tables  have  been  published  under  the  title  "Tables  of  Free  Molecule 
Flow  Functions,  Appendix  C  to  LMSC  Technical  Report  448514"  and  has  been 
given  the  number  LMSC  448514-1. 
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